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ABSTRACT 


Level Mountain is a Late Cenozoic shield volcano 
located near 131°20'W,58°31'N in the Stikine Volcanic Belt 
of northern British Columbia. The volcanic plateau with an 
average elevation of 4500',(1372m), is younger than Upper 
Miocene in age and is comprised of up to four sequences of 
alkali basalt and ankaramite flows and tuffs. The central 
region of peaks and ridges rises to a maximum elevation of 
7200',(2195m), and has a repetitive bimodal distribution of 
alkali basalt and peralkaline salic lavas and tuffs that 
Spans the period from 4.5 million years b.p. to recent 
times. Feeble basaltic vents with spatter, bombs and scoria 
apparently postdate continental glaciation. A paleomagnetic 
Study on two stratigraphic sections that span the entire 
range of volcanism samples the earth's major polarity 
reversals for the past 6 MY. 

Petrochemically the lavas belong to the sodic alkali 
basalt series and are of continental affinity. Oxygen 
isotope values near 5.6°/,.0'*O SMOW have been measured for 
unaltered rocks including basalts, salics, and ultramafic 
inclusions. This would seem to indicate that the mantle 
under northwestern British Columbia is normal with respect 
to oxygen and that the Level Mountain magmas are mantle 
derived. Corroborating evidence is available from the whole 
rock strontium isotope studies which show a total range of 
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Salic, undersaturated and ovesaturated, peralkaline and 
metaluminous rocks attests to the complexity of petrogenetic 
processes at Level Mountain. Although the two volumetrically 
important lava types appear to be primary and mantle 
derived, the presence of some rhyolites enriched in ‘80, 
granitic gneiss incluSions in basalts and tristanites and 
great variation in whole rock lead isotope ratios are 
probably all indicative of some degree of interaction 
between magmas and crustal rocks. Major and trace element 
geochemical variations and mineralogy indicate an upper 
mantle origin at shallower than 15 kbar pressure from an 
undersaturated ultramafic source that possessed an abundance 
of alkalis and incompatible elements. The mineralogy and 
chemistry of the major rock types indicate conditions of low 
oxygen fugacity and a dry gas phase composition for both 
peralkaline flows and basalts. The field relations of some 
comendite flows implies fluid behavior and seems to indicate 
a viscosity two to three orders of magnitude lower than 
would be predicted on the basis of silica content and whole 
rock chemistry. This could be partially reconciled if the 
gas phase had been halogen rich and dry. 

Calculations of the volume and energetics of the 
volcanism of the Eee Belt indicate that less than 15% of 
the modern heat flow is necessary to maintain the levels of 


volcanicity since Pliocene times. 
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CHAPTER 8. CRYSTALLIZATION AND MELTING EXPERIMENTS AT ONE 
ATMOSPHERE TOTAL PRESSURE ON SELECTED ROCKS FROM LEVEL 


MOUNTAIN 


INTRODUCTION 

A series of more than one hundred reconnaissance 
experiments was conducted to determine the crystallization 
behavior of fourteen Level Mountain lavas. Six of the lavas 
were peralkaline salics and five were basaltic. The 
remaining three were an alkaline trachyte, a tristanite and 
a phonolite. Additionally, thirty-five experiments were 
conducted to determine the solidus temperatures of six 
xenoliths of mafic and ultramafic composition. For these 
experiments, the conditions covered the temperature interval 
from 800°C to 1300°C with the principle range of oxygen 
fugacities between the FMQ and MW buffers. Twenty-nine 
experiments were additionally performed under more oxidizing 
conditions to assess the effect of variable f£0.,. For the 
most part these experiments were rock-air buffered somewhere 


between FMQO and HM. 


PREVIOUS WORK. CRYSTALLIZATION AND MELTING RELATIONSHIPS OF 
PERALKALINE ROCKS 

A number of studies have been conducted on the melting 
characteristics of peralkaline rocks. Rather different 


behavior has been demonstrated for plutonic versus volcanic 
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undersaturated peralkaline rocks (agpaites).' The most 
prominent initial observation from experimental work was 
that agpaites had much longer melting intervals than common 
rock types (Piotrowski and Edgar, 1970; Sood et al, 1970). 
The large liquidus - solidus interval has been attributed in 
part to peralkalinity and in part to high concentrations of 
uncommon or incompatible geochemical elements (Ti, Zr, REE, 
Fe**, Cl). The halogens and water in these systems may be 
preferentially retained in the melt, thus affecting an 
abnormal depression of the solidus, particularly at higher 
pressures (Anderson, 1974; Manning, 1981). Edgar and Parker 
(1974) successfully verified this hypothesis for eight 
erystallinéorocks, four plutonic ’and four tvolcanie.* The 
melting intervals were typically greater for the plutonic 
rocks than for the volcanics of similar composition. 
Presumably this was due to differences in previous cooling 
history of samples versus incorporation of halogens into 
crystals and the glass. The experiments reported by 
Piotrowski and Edgar and by Edgar and Parker were conducted 
at either one atmosphere or at one kilobar water pressure, 
unbuffered. Their two-point hydrous liquidus curves appear 
to be normal, with negative slopes (aP/dT<0). in the 
determination of the liquidus for the crystalline samples, 


' Peralkaline plutonic rocks have been termed agpaites 
(Sorenson, 1970). They are characterized by a peralkalinity 
index greater than unity (Ac+Ns in the norm), high 
concentrations’ of Pe, Ti, Zr, REE, F, Cl and unusual 
mineralogy which may include eudialite, rinkite, 
rosenbuschite. Typical occurences of these rocks are found 
at Khibina,Kola and Ilimaussaq. 
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For the four Tenerife phonolites the one-atmosphere liquidus 
temperature ranged from 1030° to 1200°C. The crystallization 
sequence at one atmosphere was feldspar followed by 
nepheline then amphibole. When pH,0 is low, clinopyroxene 
may crystallize before amphibole or in place of amphibole. 
The order of appearance of clinopyroxene on the liquidus is 
quite variable. The phonolites were of the mildly 
peralkaline type with A.I. ranging from 1.03 to 1.19. 
Melting intervals determined for three of the phonolites 
ranged from 305°C to 475°C. 

The silica oversaturated peralkaline system has been 
Studied experimentally by Bailey et al (1974). Their studies 
were anhydrous but unbuffered and conducted over a range of 
pressures from one atmosphere to approximately 1.5 kilobar. 
Their samples were pantellerites and pantelleritic trachytes 
from Bast Africa (Menghi, Longonot, Eburru). The work of 
Bailey et al (1974) was done on very peralkaline 
compositions with A.I. = 1.38 to 2.18. While the solidus 
temperature was not determined, the liquidus temperature 
showed a strong linear decrease with increasing 
peralkalinity. For all compositions studied, plagioclase was 
the silicate liquidus phase to one kilobar. Sample 83, 
Eburru Pantellerite, also had phenocryst quartz at one 
kilobar. There was something noteworthy and rather unusual 
about the liquidus curves reported in Bailey et al (1974) 
and Bailey and Cooper (1979), see Figure 8-1. The boundary 
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Figure. 8-7. 
Liguidus curves for peralkaline trachytes (KL4, KE4) and 
more Siliceous pantellerites (KE12, KW5) taken from Bailey 
et al (1974). AS presSure increases the liquidus temperature 
goes through a minimum, then increases linearly to about 
S5kbar. Depending on the ascent path, such melts could 
Gevelop considerable superheat. Note convergence of liquidus 
curves to higher P and T. 
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crystallization surfaces for alkali feldspar and quartz, are 
concave to higher temperatures with a temperature minimum at 
about 0.3 kilobar. If this is indeed the case, a peralkaline 
magma could possibly acquire some degree of superheat upon 
ascent (Waldbaum, 1971; Marsh and Kantha, 1978). Partially 
resorbed alkali feldspar and quartz phenocrysts are a common 
feature of sparsely crystalline and glassy peralkaline salic 
lavas. Comendites and pantellerites typically display flow 
morphologies characteristic of fluid behavior, which 
contrasts with the more viscous nature of calc-alkaline 
rhyolites of similar SiO, content (Shaw, 1973; Schmincke, 
1974). Both of these observations might arise from such a 
heat input implied by the liquidus determinations of Bailey 


et al (1974). 


EXPERIMENTAL TECHNIQUES 

Experiments were cinducted on the Deltech furnace under 
controlled £0, conditions achieved by the use of known . 
mixtures of CO and CO,. The samples, either pressed mineral 
powder briquettes or glasses quenched from 1300°C, were 
Sintered onto a Pt wire loop and suspended in the thermally 
uniform position of the furnace after the technique of 
Donaldson et al (1975). Experiment durations were generally 
one to three hours. Rapid quenching was achieved by severing 
the suspension wires so that the sample fell into a mercury 


filled receptacle at the base of the furnace. 
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RESULTS 

The results of the melting experiments on Level 
Mountain lavas are summarized as a series of bar graphs in 
Figure 8-2. The basaltic lavas had the crystallization 
sequence oxide, olivine, clinopyroxene, plagioclase. The 
cryStallization data is summarized in Table 8-1. An oxide 
phase crystallized above 1260°C. Compositionally these were 
magnetite - ulvospinel solid solutions. Forsterite-rich 
olivine was the first silicate phase on the liquidus between 
1250°C and 1265°C. With the exception of Ne-normative 
hawaiite, 8/25-50/6397, diopsidic clinopyroxene appeared as 
the second silicate phase. The clinopyroxene liguidus fell 
in the interval 1250° to 1165°C. Towards the lower 
temperature end of this range an intermediate plagioclase, 
labradorite to andesine composition, was on the liquidus. 
This compares favorably with the one atmosphere 
crystallization experiments conducted on alkali basalt by 
Thompson and Flower (1971) and Duke (1974). Residual phases 
such as alkali feldspar and nepheline did not appear until 
near the solidus, between 1100°C and 1050°C. Often two 
distinct forms were present in the crystallized phases. 
While the high temperature titanomagnetites tended to be 
well formed cubes or octahedra, a second dendritic or quench 
textured magnetite formed particularly in the runs quenched 
from the lower temperature range of the melting interval. 
Olivines formed as discreet prisms of equant to acicular 


habit, while clinopyroxenes crystallized as blades both 
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Figure 8-2, 
Buffered crystallization sequences and melting intervals for 
Level Mountain lavas. Comendites: (1) 9/2-99 an@ (2) 
8/25-54/6345; Pantellerites: (3) 25/5e and (4) 8/5-8/5150; 
Comenditicetmachyte: (5) 24/2 o/c and) (6) 25/7- Trachyte: 
(7) 24/hi; Tristanite: (8) LMIIIA; Phonolite: (9) LMI 201i; 
Hawaiite: (10) 8/25-50/6397; Alkali Basalts: (4) vand 31:4) 
Ankaramites: (12) and (13). T range of experiments indicated 
by length of line. Highest T occurrence of a mineral plotted 
as appropriate symbol. 
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Table 8-1] Liquidus and Solidus Data for Basic and 
Intermediate Lavas from Level Mountain 


LMI20i VocoteU" | 1050 225. | 175 245 
phonolite 


LMIIIA 
tristanite 


8/25-59/6397 
hawaiite 


PAP 
Alkali basalt 


01 
Ankaramite 


h 
ankaramite 


207 T=f ear ty | 12/0 eis 


hawaiite 


* French (1971) 
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Singly and as feathery, quench textured aggregates. 
Plagioclase crystallized both as blades and dovetail prisms. 
The crystallization sequences found experimentally 
correspond to the natural behavior inferred from the 
petography of the respective lavas; however, the actual 
crystal forms may differ due to kinetic effects (Donaldson, 
1976). The average measured silicate liquidus temperature 
for five Level Mountain basalts was 1265°+5°C, with a 
corresponding average solidus of 1083°+6°C. The average of 
the five basalt melting intervals was 182°+5°C. 

The crystallization sequence of the phonolite and 
tristanite were oxide, plagioclase, alkali feldspar, 
clinopyroxene. The oxide appeared between 1212°C and 1247°C. 
Plagioclase and andesine to oligoclase composition appeared 
above 1205°C with the silicate liquidus temperatures for 
phonolite=and  errstanite, setpatii225°>220° and \1230°225° 
accordingly. The averages for these two intermediate lavas 
are liquidus temperature: 1228°C+3°, solidus temperature: 
1045°C+5° and -mettinceentercvalts 10s (Cry -oeThe  lrqurdusrand 
solidus are marginally lower than for the basalts and the 
melting interval is the same. 

The crystallization behavior of the salic lavas is 
Summarized in Figure 8-2 with temperature data listed in 
Table 8-2. The average silicate liquidus temperature for the 
seven Level Mountain lavas was 1216° + 13°C. The average 
Salic solidus was 894°+15°C and the average melting interval 


was 308°+22°C. The largest melting range observed was for 
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Table 8-2 Liquidus and Solidus Data for Salic Lavas 
From Level Mountain 


24/1 hi WSO t10 4965. + 225 £45 
trachyte 


24/2-0/7¢ 


25/5e 
pantellerite 


25/7 


9/2-99 
comendite 


8/5-8/5150 
pantellerite 


28/25-54/6345 
comendite 


* French (1971) Plagioclase; Clinopyroxene; Olivine 
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Pantellerive 6/5 '-8/5150 (3885223 °C). The’ liquidus 
determinations are somewhat higher than previously reported 
(Bailey et al, 1979). This is presumably due to the present 
Study being anhydrous and buffered at lower f£0,. The long 
melting intervals are comparable to data on peralkaline 
phonolites and plutonics reported by Edgar and Parker 
(1974). 

Aenigmatite was the first silicate to crystallize for 
Six out of seven samples under the conditions of low f0,. 
The aenigmatite crystallized as a deep reddish brown pseudo- 
hexagonal thin platelets both singly and in "squadron 
formation", where platelets in small regions of the glass 
were all in parallel orientation, see plates 8-1 and 8-2. 
This presumably indicates regions of intermediate to long 
range order in the glass. The aenigmatite may be preceeded 
by an oxide phase, although a reaction is uSually indicated 
due to disappearance of the oxide at and below the 
temperature of aenigmatite crystallization. At the silicate 
liquidus temperature the oxide reacts with the peralkaline 
component (Ns) of the melt to form aenigmatite. The 
incongruent melting of aenigmatite has previously been 
reported (Lindsley, 1971; Lindsley and Haggerty, 1971). 
Incongruent melting has also been noted for acmite (Bailey, 
1969). Aenigmatite is followed by intermediate to potassic 
alkali feldspar and sodic hedenbergite in either order. The 
pyroxenes are prismatic while the feldspars are either 


lathlike, in the upper end of the melting range, or a matte 
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Plate 3-12 
Photomicrograph (P.P.L., 340x240microns) of a 
crystallization experiment on peralkaline trachyte 
8/5-8/5150. Aenigmatite reddish brown pseudohexagonal 
platelets occur on the liquidus for T-f£0, conditions between 
M-W and FMOQ. 
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Platese-2. 
Photomicrograph (P.P.L., 68x48microns) of a crystallization 
experiment on pantellerite 25/5e at 1201°C and log f0, = 
-8.78. Pseudohexagonal brown aenigmatite is followed at 
Slightly lower temperatures by the crystallization of an 
oxide phase. Possibly having crossed the boundary curve for 
Aent+O,=TiMt+Nds+Si0, as a result of the specific f0, 
conditions or early aenigmatite crystallization having 
lowered the peralkalinity sufficiently to shift the 
equilibria. 
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of microlites near the solidus. A bluish alkali amphibole, 
presumably riebeckite, was observed in one charge of 
pantellerite (25/5e) and in the mid-melting range of some of 
the earlier unbuffered experiments. The unbuffered 
experiments had slightly lower liquidus temperatures 
(average of five determinations 1138+35°C) with the sequence 
oxide, alkali feldspar, clinopyroxene, alkali amphibole or 
oxide, aenigmatite, alkali feldspar. In considering just the 
latter low £0, experiment, there is a great deal of 
variation in crystallization sequence and in temperature 
intervals between sequential appearance of the various 
Silicate phases. This degree of variability in the 
crystallization behavior of these seven samples would not be 
expected from the small degree of compositional variation in 


the major element contents of these lavas. 
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CHAPTER 9. THE PHYSICAL BEHAVIOR OF LEVEL MOUNTAIN LAVAS 


VISCOSITIES: ESTIMATED AND EXPERIMENTAL” 


Theoretical Viscosity Calculations 

' For lavas of known composition, viscosities have been 
predicted according to the model of Shaw (1972). These 
calculations were performed using a fortran computer 
programme supplied by J. Nicholls to give a log-linear 
expression of the form: 

ln. n(pascal-sec) = A/T°K - B 

which corresponds to newtonian behavior. This type of 
relationship has been shown to be valid above the liquidus 
temperature while below the liquidus the viscosities 
increase more rapidly due to the presence of crystals and 
the behavior is no longer newtonian (Shaw, 1969; Scarfe, 
1973; Murase and McBirney, 1973). The conversion from SI to 
cgs units, which have been conventionally used is: 

10 e exp(ln n(pascal-sec)) = n(poise) 
To characterize the expected variation in viscosity for 
Level Mountain lavas and to compare it to other British 
Columbia lavas, the theoretical anhydrous viscosity has been 
calculated at 1300°C. The range of predicted viscosities is 
tabulated as a function of rock type for Level Mountain and 


other Cordilleran lavas, table 9-1. The lowest predicted 


' Throughgout this chapter the symbol nm is used to denote 
viscosity. 
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Table 9-1 Range of Predicted Viscosities as a Function of Rock Type 


Rock Type 
Alkali Basalts and Ankaramite 
Hawaiites 


Mugearites, Benmoreites 
and Phonolites 


Trachytes and Tristanites 


Comendites and Rhyolites 


In 


0.4 to 


LeGLEO 


Aske £0 
6 69 FO 
9.0 to 


(pascal-sec)1 


2.4 
3.4 


6.4 
8.8 
10.8 


median 


127, 
(ied 


Das 
i.3 
OR/ 


1. Viscosity calculated from chemical analysis on an anhydrous 


basis at 1300°C after the method of Shaw, 1972. 
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viscosities are for ankaramites and other basalts. Melts of 
basaltic composition have a high proportion of network 
modifier cations and a less polymerized polyanionic melt 
framework (Scarfe, 1973, 1977; Scarfe et al, 1979: Mysen et 
al, 1980). The highest predicted viscosities are for 
metaluminous rhyolites with comendites being slightly lower. 
The samples from Edziza and other late Cenozoic lavas in 
British Columbia display the same range of variation as seen 
at Level Mountain. For the purposes of comparison, all of 
these viscosities were calculated using dry chemical 
analyses. The effect of water is to de-polymerize the 
Silicate melt by increasing the proportion of non-bridging 
oxygens. The calculated effect of this depolymerization for 
these relatively dry alkaline lavas is generally less than a 
few percent lowering in ln nm and the inaccuracy introduced 
by comparing lavas at a single temperature of 1300°C should 
outweigh this. The value of 1300°C was chosen to be above, 
but close to, the liquidus for all of the alkaline lavas. 
This is about 35°C above the average basalt liquidus and 
85°C above the average salic liquidus. 
Experimental Viscosity Determinations 

Three rocks were chosen to represent the basalt to 
comendite sequence at Level Mountain: hawaiite (8/25 - 
50/6397), alkali basalt (PAP), and comendite (8/25 - 
54/6345). For chemical analyses please see chapter on 
petrochemistry. Rock powders (which passed the 325 mesh 


screen) were prepared from these specimens. Experimental 
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work was performed by C.M. Scarfe at Sheffield University 
(Scarfe and Hamilton, 1980). Glasses free from bubbles and 
crystals were fused from whole rock powders, under reducing 
Gondbtions; Pinzadgas@iurnacetat 114002 =c1600°*, v¥iorCewo 
hours. Fusions were performed in alumina crucibles and the 
liquid was poured into a steel mold to make a glass slug 
whose diameter matched the internal dimensions of the outer 
cylinder of a concentric cylinder viscometer (Scarfe 1977). 
Under measurement conditions, the 8cc glass slug filled the 
Cavity between the inner and outer platinum cylinders of the 
viscometer. During a meaSurement, shear across the silicate 
liquid transmits torque to the inner cylinder which results 
in an angular displacement. The displacement is recorded as 
a light spot deflection on a linear scale. At constant 
temperature, there is a linear relationship between rotation 
speed of the outer cylinder and angular displacement of the 
inner cylinder. This corresponds to a plot of shear stress 
versus shear rate, which for newtonian fluids is a straight 
line through the origin. All experimental determinations 
were performed in air. Viscosities are considered accurate 
to #5°C intthecranges 1200A€ to 1400°C. Starting glasses for 
the hawaiite sample were prepared under both oxidizing and 
reducing conditions. The resulting difference in measured 
viscosities was very small, with the reduced aliquot being 
slightly less viscous at any temperature. This could be 
related to enhancing the behavior of Fe as a network 


modifier under reducing conditions (Fe’** rather than Fe**). 
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Results for the three Level Mountain lavas are presented in 
Figure 9-1. The apparatus was calibrated using NBS Standard 
710, a soda-lime silica glass, for which the 
viscosity-temperature relationship is accurately known. 
Curves through experimental points are fit by hand. Dashed 
lines correspond to theoretically predicted viscosities. 
Viscosities calculated using the method of Shaw (1972) agree 
with experimental determinations to within 0.25 log units, 
with the predicted values being lower than measured values 
in all cases. A linear relationship exists between the 
logarithm of the viscosity and the reciprocal absolute 
temperature in the range 1300° to 1400°C. The temperature 
dependence of the viscosity can be described by an Arrhenius 
relationship of the form n = n,e exp(EnN/RT) where n, is 
constant for a given melt, and En is the activation energy 
for viscous flow. Activation energies for the Level Mountain 
lavas are in the range 98 kcal/mol for the comendite to 54 
kcal/mol for the hawaiite. At any temperature above the 
liquidus the viscosity of the comendite is about three 
orders of magnitude greater than the basalts. For example, 
at 1300°C the comendite has 1n n = 5.3 poise as compared 
with Yin Wp = 2. 050and 12.160 for the basalts: (in 'n = 924, 3.6; 
2.3 pascal - sec respectively). Despite the fact that 
theoretical and experimental data agree in predicting a high 
viscosity for comendite, this does not account for the 
apparent fluid nature of comendite flows at some locations 


on Level Mountain. Both experimental and theoretical data 
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Figure 9-1. 
Plot of reciprocal absolute temperature versus log viscosity 
(in poise) comparing experimentally determined viscosities 
to Shaw model predictions, taken from Scarfe and Hamilton 
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relate to a degassed dry sample. Under natural conditions of 
emplacement and extrusion, dissolved volatiles in the 
comendite melt, particulartily fluorine and chlorine (Bailey 
and MacDonald, 1975; Scarfe, 1977), may have reduced 
viscosities considerably below 10° poise. Degassing could 
eventually increase viscosity to the point where flow will 
eedset (Nash ~and Bvatish!1977)..Phesether(possibilitysis that 
comendite may acquire some degree of superheat (Marsh and 
Kantha, 1978) due to a peculiarity of the liquidus geometry 
in PT space and ascent path (Bailey et al, 1974). The only 
problem is that for the viscosity-temperature relationships 
presented here, an extrusion temperature of 1500°C would 
Still give a viscosity too high to explain fluid features, 
lava tubes, thin flows (less than 3m thick), long low volume 
flows, lineations, flow folds and thin dykes (less than 1m 
thick) which have been observed in the field. Comendites 
Should have a more polymerized melt than basalts due to high 
silica contents; “highsR=(Sit+Al) /O: (atomic) (Shaw; 91965? 
1972) and low ratios of non bridging oxygens to 
tetrahedrally coordinated cations, principally Si and Al 
(NBO/T, Mysen et al, 1980). The precise effects of high 
Fe/(Si + Al), (Na + K), (Na/K) halogen content, and eruptive 
temperatures, as compared with these values in calc-alkaline 
rhyolites of similar silica contents, are unknown; however, 
all of these characteristics could contribute to a lower 
viscosity (Schmincke, 1974). Also there is probably a 


pressure effect which contributes to decreasing the 
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viscosity of salic melts due to a collapse of their 
one-atmosphere 3D polymer structure at higher pressures 
(Virgo et al, 1979). Characterizing the viscosities for the 
Level Mountain lavas is important to an understanding of 


their flow properties, extrusion rates and nodule transport. 


CONSTRAINTS ON THE RISE TIME OF LEVEL MOUNTAIN LAVAS 


Introduction 

This section is devoted to placing physical and fluid 
mechanical constraints on the rise times and emplacement 
velocities of the various lava types found in the alkaline 
association at Level Mountain. Settling calculations have 
been performed for selected flow-nodule pairs. Some diapiric 
rise and heat transfer calculations are also presented. The 
data used in this section comes from a variety of places. 
Liquidus and solidus determinations were made on glassy 
charges suspended in a Deltech vertical tube quench furnace 
equipped for gas mixing. Oxygen fugacity was controlled to 
lie between FMQ and MW buffers using mixtures of CO and CO, 
after Dienes et al (1974). Viscosities are calculated at the 
liquidus and, when available, experimental values are used. 
Due to the dry nature of these lavas (See chapter on 
petrochemistry), the effect of wet versus dry viscosities is 
small, especially when compared to other assumptions made in 
model calculations. All of the liquidus data, eleven values 


covering the entire range of lava compositions, were plotted 
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as 10‘/T°K versus 1n n(SI units) to give the following 
linear relationship at the liquidus: 

Inen l= eraons Ghet GEo Ka bk7 90 
with a coefficient of determination of r*= 0.79 as 
determined by least squares linear regression. This 
empirical expression was used along with the Shaw equation 

iia bee AA TAK AaB 
for each lava to predict liquidus temperature and viscosity 
where such data wasS not otherwise available. Reasonable 
liquidus estimates were also obtained using this expression 
and other B.C. alkaline lavas, the McKinney basalt, and some 
calc-alkaline andesites as well. It is conceivable that an 
expression relating viscosity and temperature should exist 
for a given volcano. This could be due to the mechanism of 
magma genesis for a related suite or perhaps due to the 
effect of common plumbing. Estimates of liquidus temperature 
and viscosity could also be made by comparison with other 
lavas without recourse to speculative and empirical 
expressions. The determination of the significance of this 
type of expression would require further research with 
additional data on viscosity and liquidus temperatures. The 
rational for using such empirical ‘calculations on Level 
Mountain lavas was to be able to use values for actual host 
lava-nodule pairs. 

Xenolith densities were analytically determined uSing a 
weight and volume displacement technique. Lava densities 


were calculated at the liquidus temperature after the method 
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of Bottinga and Weill (1970). The Bingham model of Sparks, 
Pinkerton, and MacDonald (1977) was calculated using their 
empirical graphs and SI units. Other calculations were done 
for convenience in cgs units. 

Nodule Transport and the Ascent of Magmas 

The majority of lava flows on Level Mountain are 
aphyric to fine grained and relatively free of phenocrysts. 
Newtonian flow models should be adequate for these lavas, 
however Sparks et al(1977) made a case for phenocryst laden 
alkaline lavas behaving as a Bingham fluid. It is 
conceivable that up to several flow volumes of alkali magma 
could be transported through the crust as crustal rich 
Giapirs. Information relating to the physical properties of 
lavas used in model ascent calculations is presented in 
Table 9=2'. 

Petrographic and field relations suggest that most 
lavas were at or near their liquidus temperature at the time 
of eruption. Pleistocene basalts from Meszah Peak contain 
xenoliths of anorthositic gneiss, troctolite and gabbro. Of 
these the troctolite and gabbro are considerably denser than 
their host flow. Similar basalts of the adjacent Stikine 
canyon occasionally contain lherzolite nodules (Littlejohn 
and Greenwood, 1974). At three locations on Level Mountain 
peralkaline trachytes were found to contain armoured nodules 
of ultramafic composition and a poSitive density contrast 
with their host flows. A two meter wide tristanite dyke 


contained 5 cm long gneiss fragments with a positive density 
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contrast. Further description of these nodule flow pairs is 
presented in Table 9-3 and in the chapter on petrography. 

The existence of xenolith fragments in certain flows 
places two constraints on the rise times of these magmas. 
Where the nodule density is greater than the host lava and 
the nodule is sufficiently large, a slow ascent velocity 
would permit the nodule to settle out. If the solidus 
temperature of the nodule is below the liquidus temperature 
of the host magma, aS it appears to be for most examples 
given here, the nodule should react, melt and possibly 
disaggregate. The rates of these phenomena depend on nodule 
Size, relative thermal conductivities, grain size, fabric, 
chemical diffusivities and heat capacities of nodule and 
melt, as well as on the melt viscosity and flow regime. For 
the nodules to arrive intact at the surface, their residence 
time in the melt must be short enough to prevent either 
settling out or destruction. 

For flow models involving a viscous compressible fluid 
and low concentration of spherical particles, the settling 
behavior changes as a function of frictional drag and 
Reynolds number. The Reynolds number is given by Re=(ZrVe/n) 
(Bird, Stewart’ and Lightfoot, 1962). Bor Re <<h.0 the 
frictional drag is effectively a linear asympotic function 
of Re and the Stoke's Settling Law behavior is valid. 

For intermediate Re (1.0 to 1000.0) the relationship is a 
power law, where the friction factor (f) is given by f= 


1.85Re exp(-0.6), reflecting the existence of a turbulent 
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flow around the equator and in the lee of the particle. For 
high Re the friction is nearly independent f = 0.44 and the 
type of settling is in keeping with Newton's Law- The cutoff 
of Re = 2000 marks the onset of turbulent flow. Particle 
Reynold's numbers for relatively aphyric melts are typically 
smaller than 107° for salic compositions and between 2x107' 
and 2x10‘ for basalts. For these crystal and nodule poor 
cases Stoke's law settling behavior should be reasonably 
valid. The greatest uncertainty is not the choice of the 
settling law, but the assumptions of viscosity and density 
contrast. For multiphase flow of crystal and xenolith 
mushes, the disperse particle theory no longer applies and 
particle boundary layers and particle interactive forces 
come into play. An alternate way of stating this is that in 
phenocryst laden melts, there is significant crystal-crystal 
interaction which impedes flow. Bingham model ascent rates 
have been calculated for the Level Mountain examples 
according to the models of Sparks et al (1977). For the 
basaltic magmas a rapid ascent time is still implied. The 
Bingham model is not very informative when applied to the 
Salic magmas. The minimum ascent velocities for nodule 
transport, rather than settling, imply maximum rise times on 
the order of tens of thousands of years. Upon consideration 
of cooling and crystallization processes and thermal 
conductivities of the materials in question, this is 
sufficiently long to cool salic magams by conduction alone. 


If ascent were this slow, plutonism rather then volcanism 
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would have occurred. 

The model of Takeuchi et al (1972) assumes buoyant rise 
of melts through a lithosphere of greater density. In their 
model the minimum velocity for magma rise is given by 

VPHApgda-/ 12M, 
where d is the characteristic dimension of the conduit. 
These calculations have also been presented in Table 9-4 
uSing an hypothetical conduit dimension of 10m, which is 
typical of high level vents on Level Mountain, and an 
assumed crustal density of 2.9g/cm*, after the crustal model 
for the Intermontane Belt in Southern B.C. (Keen and 
Hyndman, 1979). Again the assumptions regarding geometry and 
viscosity have far greater effects than actual magnitude of 
density contrast. In all cases this buoyant magma rise model 
predicts velocities that are orders of magnitude faster than 
the minimum velocity limits set by nodule transport. This 
model presupposes a magma-filled conduit extending through 
the crust. The driving force is the magmastatic head 
(difference in lithostatic and magmastatic pressure at 
depth). For basalts this magmastatic head may be as high as 
a few kilobars (Maaloe, 1973). The model does not account 
for frictional effects or the pressure required to hold the 
conduit open. All limitations considered, this model should 
put a maximum constraint on ascent velocity. In this case 
the velocities calculated for the hawaiite and the "typical 
basalt" are 20 to 30% of the expected seismic velocities 


(Vp) in the crust. Velocities of this magnitude are also 
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expected from the arguments for magma emplacement in the 
crust governed by fracture propagation (Anderson and Grew, 
1977). Using the model with its buoyancy drive mechanism, 
all types of magmas encountered at Level Mountain can rise 
from mantle depths in a matter of hours. 

Rise times can also be estimated by calculating the 
time required to heat a nodule through its melting point. 
This is essentially a "Newton's Law of Cooling" calculation. 
Several assumptions and preliminary observations should be 
noted. In Clark (1966) it is shown that refractory minerals 
have higher thermal conductivities than melts of higher 
Silica content. In keeping with this observation, an 
inclusion made up of refractory minerals will tend to have a 
very flat thermal profile as it heats up in a melt. Forced 
convection of the magma and turbulent flow will tend to 
transport heat to the nodule faster than assumed in these 
calculations, which will diminish the residence time to 
melting. For any positive density contrast (nodule greater 
than melt) some crossflow must exist so that this time 
estimate will be a maximum. Spherical nodules have the 
highest volume to surface area ratio so that any other 
geometry will enhance heat transfer and diminish time to 
melting. For this model the characteristic time is the 
duration of the heating process, only up until the nodule 
solidus is achieved. For real nodules once the solidus has 
been reached there will still be a finite time of melting 


associated with kinetic effects such as diffusion and 
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reactive melting. In this model, for steady state 
conduction, the rate of change of temperature with time is 
given by Tuma (1976) as 
dT/dt = -hcaTA/VpCp 
where he = heat conduction coefficient, cal/cm?-sec-°C 
AT = initial temperature differential, °C melt - 

nodule 

A,V = surface area and volume of nodule, cm’, cm? 


respectively 


9 density of nodule 


iH) 


Cp heat capacity of nodule 

The heat Ses density and geometry of the nodule 
are assumed to be constant up until melting. To extend the 
model through melting, the heat capacity would have to be 
replaced by a term including the heat of melting and any 
excess heats of mixing associated with diffusion processes. 
For a spherical nodule geometry A/V = 3/r. Newton's Law of 
Cooling states 

AT (tL) = AT. eeexp(aat) 

where a = (dT/dt)/AT.. 

The only remaining quantity to be evaluated is hc. Kreith 
(1973) gives he = kN/D where k is thermal conductivity and N 
is Nusselt's number, which for Reynold's numbers Re <10.0 is 
approximately equal to 2.0. Values of k for minerals and 
some rocks are published in Clark (1966). A sample 


calculation is presented for a general alkali basalt and a 


lherzolite nodule assuming the basalt to be at a liquidus 
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temperatunevtofui2257°C;)itherlherzolite initially at<1045°2, 


withaa Bsolrdustakh12452@4(COaA19. tkhbaryavYoder,+ 1976) .4Then: 


het=aO sO t9i2 415 =10s00b47ecal em? ssecn°C 
a€T/dtgam(<0500147)3(210)/755( 943 56940) 
=6.48x10°* °C/sec 
C=O UOK Um erie) =9 (Occand: Te =. 240i, 


Finally the heating expression can be rewritten in terms of 
VIto gives 

Pee ne GONE t 6.) / a 
which with substitution of the appropriate values, 


t 


InCO0V2 10) /—3.08x 10° * ‘gives 


ie 


Jo xe tO Seconds, or 11.43 days. 

Assuming the lherzolite to have become incorporated at a 
depth of 30km yields a minimum ascent velocity of 3.039 
cm/sec. Depending on the assumed depths of incorporation for 
the inclusions the model ascent velocities are generally 
between those estimated by Stoke's Settling and buoyant 
rise, see table 9-4, 

Carmichael et al (1977) present an ascent model whereby 
Spherical magmatic diapirs rise buoyantly, yet frictionally 
retarded, through a very viscous upper mantle. The 
relationship they derive is R*? = 3Vn/gAp, where R is the 
effective radius of the magma diapir, V is the ascent 
velocity and nm is the viscosity of the upper mantle. The 
value of n, derived from isostasy, and typically used in 
geodynamic arguments is 107' poise (Cathles, 1975). A few 


trial calculations, using reasonable density contrasts and 
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assuming a value for either V or R to solve for the other, 
reveals unreasonably large radii, of hundreds to thousands 
of km, or unreasonably slow velocities, which for depths of 
30 to 100km imply rise times of tens of millions to billions 
of years. The Cordilleran lithosphere is thinner and the 
upper mantle is softer, with lateral variations in the 
effective viscosity spanning two orders of magnitude at the 
100km depth range (Ranalli, 1980). The diapiric transport 
calculation can be rearranged to estimate a viscosity for 
the upper mantle that would be volcanologically more 
reasonable. Consider a peralkaline trachyte flow, PBE from 
Meszah Peak. The volume of the flow is 4.8x10°~*km* for an 
equivalent spherical radius of 105m. Assuming a slow ascent 
velocity of 1.16x10°* cm/sec, and densities for mantle and 
melt of 3.3g/cc and 2.4g/cc respectively, results ina 
calculated estimate for n(mantle) of 2.78x10'* poise. For a 
hypothetical diapir volume of ikm*, a velocity of i10°-* 
cm/sec and the same density contrast, the value of is 

1013x 10? iporse. aThrs begins to hbookeplausible forlsethe 
Situation of a partially molten upper mantle and transport 
of diapirs having about 10 flow volumes worth of magma. For 
the ascent velocity of about 10° ‘cm/sec, it would require 
approximately 42,000 years to traverse 13km of upper mantle. 
The evidence from age dating attempts suggests that the 
major eruptive periods are separated in time by tens of 
thousands of years. This coincidence lends support to the 


buoyant diapir concept and calculations discussed. 
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Assuming that all of the Level Mountain lavas rose as a 
Single magma diapir from a hypothetical depth of melting and 
segregation of i105km to 30km. The total volume of volcanics 
is 1500km* for an equivalent diapir radius of 7.1ikm. As the 
lavas are predominantly basic, a more reasonable density 
contrast is 0.6g/cm*®. For a normal upper mantle viscosity 
(10?'poise), the calculated ascent velocity is 9.89 x 
10-*cm/sec and the traverse time to the base of the crust is 
7.58 x 10'?sec or 2.4MY. The age for the lowest flows 
Sampled in the Little Tahltan Canyon is interpreted from 
paleomagnetic evidence to be about 6.4MY. The combination of 
these two ages, 8.8MY, should be close to the origin time 
for the Level Mountain magmas. In Hamilton, Baadsgaard, and 
Scarfe (1978) a whole rock Rb-Sr age of 9.4+1.0MY was 
reported and interpreted as possibly representing some 
precursor thermal event in the upper mantle. The coincidence 
of the reported and the modelled date may lend some Support 
to the mantle diapir model. 

Carmichael et al (1977) in using the same model, and 
reasonable model parameters for alkali basalts and 
tholeiites, derived an apparent viscosity for the upper 
mantle of 5.3x10° poise. They concluded that, while this is 
unlikely to represent the viscosity of the whole upper 
mantle, it could represent the viscosity of a thin 
lubricating sheath of partially fused mantle material 
Surrounding the diapir. They also did a conductive heat loss 


calculation, inferring a drop of about 2°K for flow sized 
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bodies ascending at typical nodule settling velocities. 

Carmichael et al (1977) also calculated convective heat 
loss for these magma diapirs with thin boundary layers. For 
diapirs in the volume range 0.03 to 0.185km?, the boundary 
layers were all less than 40cm thick, which is one per mil 
of the diapir radius. The temperature losses associated with 
laminar boundary layer convection were 21° and 7° 
respectively. The losses corresponding to a turbulent 
boundary layer were higher, 86° and 39° respectively. Using 
the appropriate physical parameters for the Level Mountain 
hawaiite flow, 8/25-50/6397, and a model transport from 
100km, the laminar convective boundary layer model gives a 
temperature drop of 3° while the turbulent case was 23°C, 
Convective heat losses offset by proportional 
crystallization could be a significant driving force and 
mechanism for differentiation even at mantle depths. This 
might be a means to derive large volumes of hawaiite from a 
more primary composition while still in the mantle depth 
range. 

In addition to conductive and convective heat losses 
some other types of heat changes need to be considered. The 
adiabatic gradient is given by Yoder (1976) to be between 
ORAtaiid OC, SAG km AITRESeis lequiva Vent Bed 207 *t oe sOe. for a 
100km path. For isenthalpic decompression, Waldbaum (1971) 
cites a value of adiabatic heating on decompression which 
may be as great as 22°C per kilobar. Clearly this is a 


mechanism to offset convective heat losses and even perhaps 
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to provide superheat for an ascending magma body. This 
mechanism could provide sufficient heat for a magma to 
asSimilate significant quantities of wall rock, whether in 
the upper mantle or crust. Here the ascent velocity would be 
the crucial parameter. Again it is pointed out that any 
rising magma must undergo decompression and, if this can 
generate heat and drive assimilation processes, it may be 
difficult indeed to transport and erupt a primary magma. 
Thermal changes due to viscous heating and transportive 
friction may also be significant but their effect is thought 
to be less than adiabatic and convective effects (Carmichael 


et al eeas19 Tu. 


VOLUME ESTIMATES OF LEVEL MOUNTAIN VOLCANISM 

Volume estimates have been made for thirty-five typical 
eruptive events from Level Mountain. Estimates have been 
made for single eruption flows, tuffs, domes, and dykes 
representing the entire range of compositional types. 
Although estimates cover the time span of volcanic activity, 
the younger events are more heavily represented due to 
preservation and mappability. Dyke and vent pipe thicknesses 
were estimated at the base of the plateau as this was the 
volume of material added to the volcanic pile during an 
eruptive event. Volume estimates are presented in Table 9-5 
along with characterizing information for locality, map unit 
= age, and composition. The entire range of volumes is 


Gallen 10" om? SC0.6 12k) to. 02x10" m2 ecoeeOtkm?:). Typical 


a aaa wt rate soi at J Tess 


PEEVE IDG: : ie im, rent aint bis Tavis say 
ae | cai aetey tne day re wiya LE ‘fags heg sale 


, 6, pe Fale sites raaater’ afd, rag Hebi Ceivesom Fhe 
ae 


ory n ] 


20h Weed ‘gene om. an thsesain ae stily aa 
MOA ANG. See Bhat Ss aeeM ABE ong | fluo: a. 


Ns; 


atone peas LS! he) a f RANE: Tr hE Ce Te Py 


roahay: Miao ig oak Pear 20 slokom® 
isa? JUG) Gear wu ke 32 aR aaa oe 


tt $bRe no leesaaneseS Oe TSebaG eee ion 


232892 es med TELM LES 76 Oox: eat 


Rem Yreme ius ‘Soa Iyvoyearesd on Sbdane er 
1 at 
‘ {gtaAone Ss ey “S26 ha 43 ¥ 3 sub 2 
~ = i’ —FrPY ry 2 rt — ‘ wg hi a ad ioala Yom: 


hr py t- 
Me lwac a0 nae lavas: ne aaveamtinad a 
r=—y tel) ceeeet Saba npeey sy.an co oe van oh to 


a4 ottsay eLeseuan| ie mos? aznaee, 
ae 


Opa ashen Rae 5. euaed. mG ‘a ial nite si as 


t + 


as 


1S Patan Sen igen BLOM SIN Po Sagi 


1 A alias gis ‘wang? ody ae sig: ea 48 ba heme’ 


; ¥ eet 
sAgeD nv aaidnsea i, ‘eh jeeasians ‘egboy ~atore. “ 
5 aie 


pi, pies serena Hh 
ia ines oes +) bain: 


, 


ns. 
vie bry ee 
7 i ry : 
a. (ma hoa? Ohe 


. a a 


358 


LOLXB8L°¥8S0'} 
S33NY B SmMOLS DILES (614) BAY 
o1;seg 


.W4% ZOO’ 800° 
SittS *S55NY ‘SMOLS (8) Bay 


«U4 89° OFFI 
se7Ap ‘849038 ‘sSaucd (g) Bay 


cW4 pOO' F710’ O 


(€4) Bay 
SMO, 4 


.W4 7B OFbKTH'O 
(9) Bay sexAQ 


900° 


89E° 
O10" 
cEO° 


cEO° 
E60" 
8cS" 


SECO CIOS C CS 


cOXLE IF 
+O1X08 bp 
101x08'p 

yoee 
OlXTb 6 


vOLXEE GS 


sObXZ6° 


- 


sOLXtG 4 
»ObX80'8 
«04X00 ' 
»Ob}XO8' 
sObXbo th 
sObXt1 9 
sOLX40'% 
»Ob X09 4 


104 XO" 


— 


»ObX88°E 
sObXvG + 


20X79 L 


:ObXbb'S 


10X07" 


- 


eOlXZO'4 


»Ob XO9° 
eObLX9L” 


7 


sOlX19°Z 
eOLXxGZ'9 


+OlLXB9'€ 
»ObX65'6 
.ObxXGI'€ 


sOLXLE'E 
1ObXG2'6 
vObXB8Z'S 


(64) cW SWNLOA 


SUGAR 


Le@ssAqedAH pue eA}3dN4Q ULB UNOW 


Aat ren eAnyonyey 
eps AG NyUOPIPNG 


te /Ay eeEg 


48eq YeZSow 

gO 383 e6piy 

4p Ay eeqg AYO 'N 
“Eten A@LLeA “S 
4eeud 

Ay eeqg 4O HUO4 
eAnyonyey 
/AX+e04 EPLALA 
eAnyonyey 
/A¥¥B09 EPILALG 
4#@G@dg YURZSeW 4O 

“SS A@Bi Lea Bu, Buey 


uyeyuDd 1 BuzUEeD 
"Bed 

uezsew 4O ‘N 

2 pUWNS 

1e@d YRZSOW 
uoAue) 

UBILYeL B143h7 
u;6uew 

Nee{Yeld U4e}semM 
Aap ten eAnyonyey 
e6p+y 4Yy4ON 
NBOBlg U4 Seq 
e6pi4y 

NZUOPIPNG 4YANOS 
u;6uew 

NBOBLd U4eIsSaM 


uoAueod 

UBILYRL B12317 
up 6uew 

NeeIBl dg U4aISemM 
488g YeZSeW JO 
“Ss A@atten Bu; Buey 
NWOLeL dg Use ,sSeW 
(MS) 8487 
WNYDJOyH 814347 
uoAued 

UBILUBL O14347 
(mS ) 

NEGEV id LLeYPB4M 
NeBBid 11e4Pe4m 
NBBBid LLBUIBIM 


(aN) 

NBBISld 1 LeYYeuM 
“yO Ay}eeG AYO4 M 
NB@BIBid “3°N 


Ay} (e907 


090S/76 
Lt+S/uy}1oo0e4 


quog OZIW7 sewoqg 


2+409S 8 $3N1 G 


00L9/6L 

B}4OOS @ 45NL 
67-11/8 

Pray 
St6p/1-01/8 
331 

OGE9/09 

ewog 
0689/7S-S2/8 
se7Ag Asset 
38d 

mold 

VIIIW)’ vy-91/8 
ang 
08eS/7S-Sz/8 
4903S 

18d 

pepudog B @yuoos 
dVd 

MO({4 YWeUuWN|OD 


AOL4 “weuWN,oD 
S1€9/08 

MO, 4 4euWNL Od 
MOL4 PEPpudod 
t/ve 

MO, 4 seuWN,oD 
4d/8/GL 


Ovce/9 
MOld VV 


mold VV 


€819-8119/69 
}U@A B},4uodS 
eANL 


40 eAny 
SMOL4 
SMOL4 3 SexAQ 


MOL4 @dpd Ua 
eayAg 


61 /9t 


Ples4 BAB uO 
eyAg quwoos 209 
64AQ 
8AQ 


eunyeay 


eyAYyoReuy “9 


eyAudeu, *9 


834 1OAua 

@}} 48, 1ezUeg 
G8y+4eL1e Ue, 
ey yAYdeu, 

8}; PUeWOD 

ea yAyoRus, ‘gd 
8} +}UR Spy 
ey AYuoeu, 

By} }eMeH 
+(eseq 

84) }8MeH 


Oy}, emMeH 
+1 eseg 


F1SSbq FLEHLV 
@,990ug 

33N1 82,U06e eg 
#185889 FLENXLV 
8}; weueyHUY 


Oy} peMeH 
+1eseg 


+18Seg FLBALV 
@} }weueyuy 
2} ;weueyuy 


8} }weueyuy 
@} pwwBuBeyHUYy 


8} }wWeueyxUy 
F1BSSG PLB 
4leseg FIBNV 


adA, eubew 


eg 


eL 


v 


+4UN dew 


L®Ae7] ewoS 4O SB}BwW}YSZ BWNLOA “G-G EL Qe, 


eee 


tanto Tupcgys 
‘al? Nie ; 


6¢7/ 4 


‘Sel 
*Qy he 


tS) Wao 


‘aad 4 


rb ois“) eneerrn 


a | gies 


une ait fe 
ye Heras wy: o)1) 
‘ (t/a 


Lv) ioe. ea wy 
ia ures 


‘2 ised) Mas vi? 
we 


—_ Ts ..- 
oe 


> ‘Cam e 


e044 terme o9 
. 48 


* Jethro a) “—< 
- fe 


26 Oy yeeata 
(Os Ke te 


, oro 

iy 

" rar 

ecae\? a de 
J 


wei | f Trae} 

’ hj ’ +3 
‘ Yt Ta a eek 

« Tir3\* 
“ La ‘ ‘4 i « 

i wt 3m Wye ¢ 
ia ty Or a 
Garth 


el tretpe=t 
arenant 4 
orga 

} pene 
tie iD 
sare et 


uf i (eet 


ees 


Suagae 


a. 


ey - 


Se 


359 


historic basalt flow volumes from Hawaii and elsewhere are 
given as 0.03 and 0.185km* (Carmichael et al, 1977). The 
average volume of five post-glacial basalt, basanite and 
nephelinite cinder cones from the Cariboo, Chilicotin and 
Stikine regions of B.C. is 0.0134 km*. A large Pleistocene 
basaltic tuff cone near Jacques Lake (Campbell, 1961) has a 
volume of 0.126km*. Compared to these, the Level Mountain 
events appear to be of ordinary size. The volume estimates 
have been subdivided into basic and salic types as well as 
into flows versus hypabyssal intrusions. It is apparent that 
at Level Mountain the dykes, stocks and domes are an order 
of magnitude larger in volume than the flows and ejecta. 
Flow volumes are taken to be single eruptive events, while 
dykes and stocks may be intermediate in size between 
eruptive events and eruptive episodes. There is no 
Statistical difference in size for basic versus salic types 
in either the eruptive or hypabyssal group. This would seem 
to indicate that eruptive volumes are determined by common 
energetics of magma rise or by a common set of volcanic 
plumbing, rather than by the relative volumes of melting or 
differentiation events. 

The volume of eruptive episodes, taken to be synonymous 
with the sequential map units, have been estimated using 
topographic maps at 1:50,000 scale for elevation control and 
a grid overlay for area. Volumes have been reconstructed to 
include materials eroded form the plateau margins and 


youthful alpine glacial valleys. Model densities have been 
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assigned to each unit on the basis of average physical and 
chemical properties for rock types present. Time spans for 
map units have been estimated from all available sources of 
field and laboratory data. The average time intervals 
between flows have been calculated for a typical flow volume 
of 0.01km* in size. This data is Summarized in Table 9-6. 

There have been approximately 85,000 flow sized 
eruptions at Level Mountain over a time period of about 6 to 
7MY. The lifespan of the volcano is clearly divided into 39% 
initial plateau building and a 61% later stratocone 
development. The average time span calculated for flows of 
the plateau building period is 43 years while for the 
statocone development phase the interval is 226 years. Field 
mapping and geomorphology indicate that the activity was 
likely to have been more intense with widespread 
contemporaneous eruptions separated by periods of hundreds 
of years of quiescence. The figure of 43 years for the 
plateau period seems reasonable as there was insufficient 
time between sequential flows for extensive soil development 
Or vegetation to have occurred. 

Of the total 857.6km° of volcanics, 80.8% was basaltic 
while 19.2% was salic and dominantly peralkaline. The 
plateau building site was wholly basaltic. The stratocone 
Stage, with a repeated bimodal process, was 85% salic in 
character. 

In addition to the Level Mountain Volcanism, the 


volumes have been estimated for the other major centres in 
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the Stikine Volcanic Belt. This information is presented in 
Table 9-7 along with some estimates of dimensions pertaining 
to the source region. The total 1710km*® of volcanics 
outcrops over an area of 8629km? in the Stikine Volcanic 
Belt which has a total area of 70,000km°. The volumes and 
areas of the eleven major volcanic centres of the Stikine 
show a logarithmic distribution. Level Mountain and Mount 
Edziza are the largest in both area and volume. Using the 
total volume estimated and a time span of 5MY, it is 
apparent that the Stikine Volcanism is small in comparison 
to the Columbia River Plateau or the East African Rift and 
that the extrusion rates are lower by two orders of 
magnitude. 

The volume data on the Stikine volcanism can be used to 
evaluate two possible source models. The two models are (i) 
partial melting of a slab of upper mantle and (ii) upper 
mantle diapirs below each of the eleven major centres. The 
slab hypothesis (using the surface outcrop area as an 
estimate of source region extent) implies the removal of a 
layer up to 180m thick beneath the volcanic centres. If the 
volume is derived over the extent of the Stikine Volcanic 
Belt, the removed layer may be as*thin as 23m. When these 
Slab calculations are performed for each centre, the 
calculated layer thickness varies from 76m to 532m with a 
mean of 235m and a standard deviation that is 33% of the 
range. Assuming that 5% to 25% partial melting accounts for 


the volcanics, and that this volume may be efficiently 
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removed, implies a thickness for the melting region of 
between 0.94km and 4.70km for the mean figure quoted, or 
0.092 to 0.460km for the minimum case of partial melting 
below the entire belt. When the eleven estimates are 
considered as a population, the distribution is normal and 
85% of the entire Stikine Volcanism is accounted for within 
a thickness variation of +0.5 standard deviation of the 
mean. All of the thickness estimates lie within +2 standard 
deviations of the mean. From these model data it appears 
geologically reasonable that the Stikine Volcanism could 
Originate from a thin zone of partial melting. 

For the thermal model presented in this chapter, the 
partial melting zone could be as thick as 6km. USing this as 
a maximum estimate, and keeping the same levels of partial 
melting, an estimate of melt extraction efficiency can be 
made. A relatively inefficient melt extraction process has 
been inferred from studies of lherzolite textures (Waff and 
Holdren, 1980). For the case that melt generation is 
localized, the efficiency of melt extraction is between 
78.3% and 15.7%, while for the case. of a general partial 
melting zone beneath the entire Stikine, the melt extraction 
efficiency is between 7.7% and 1.5%. In either case a period 
of partial melting results in a considerable amount of 
material that may be left behind. For basaltic partial 
melts, the melt that stays behind could generate regions of 
gabbro in a more ordinary peridotite upper mantle. The idea 


of inefficient melt extraction is important to the concept 
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of upper mantle differentiation but has been proposed on 
other terms (Mysen and Boettcher, 1975; Lloyd and Bailey, 
1975; Boettcher and O'Neill, 1980). Sufficient 
preconditioning of the upper mantle to generate gabbroic 
pods (plagioclase plus pyroxene) could result in the 
derivation of more salic magmas in subsequent partial 
melting events. It is noteworthy that the three largest 
centres of the Stikine have both the long spans of activity 
and the salic lavas one would expect from such a two stage 
process. 

The second source model for the Stikine postulates 
separate diapirs below each of the centres which, 
SuccesSively tapped, could account for the cumulative 
individual volumes. For a spherical diapir model, the radii 
so calculated for the 11 Stikine centres are normally 
distributed. The mean radius is 2.61km and the range is 
1.06km to 5.90km with a standard deviation that is 31% of 
the range. Unlike the slab calculation, only 12% by volume 
of the Stikine Volcanism is accounted for by average sized 
diapirs. Another way of stating this conclusion is that, if 
diapirism explains the Stikine volcanism, the diapir size is 
highly variable. This implies that 88% of the volcanism must 
be accounted for by two mega diapirs of extraordinary size. 
Therefore the first source mode (slab partial melting) in 
localized or widespread regions seems to be the more 


reasonable hypothesis. 
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THE ENERGETICS OF MELTING AND RISE 


The Energetics of Volcanism at Level Mountain and for the 
Stikine 

Volcanism can be considered as an energy dissipation 
process. The thermal requirements to initiate and maintain 
this process include mantle preheating, phase change energy 
associated with melting, and the work performed in magma 
transport to the surface. The total energy released at the 
Surface can be considered as a contribution to the total 
heat flow. 

AS a monitor of the Level Mountain volcanism, the total 
energy dissipated will be considered as a function of time. 
The timing comes from magnetostratigraphy and limited age 
dating control. Mass eStimates for successive units were 
derived from field mapping and laboratory density 
measurements. Three energy terms are considered. The cooling 
of lavas from their liquidus temperatures to surface 
conditions is expressed as a heat capacity intergral, using 
heat capacities for typical flows from each unit as 
calculated in ROCK. Due to the small variation (+ 50°) in 
liquidus temperatures, regardless of lava composition, the 
value of 1200°C was used for all calculations. The heat of 
fusion values apply to liquidus temperatures for specific 
compositions considered representative of each unit. The 
transport work is expressed as a potential energy term (mgh) 


for the total mass of each map unit and the constant depth 
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of 35km. These data and calculations are summarized in Table 
9-8. The largest energy contribution for each unit is the 
cooling term which is about five times the heat of fusion. 
Map units 5a, 6a, and 7a-c are dominated by salic lavas. The 
energy components of salic lavas are distinctive. Due to the 
low heats of fusion for salic lavas, the lift energy 
required for magma transport exceeds the phase change 
energy. For the remaining units, dominated by basaltic 
compositions, the reverse is true. 

A graphical representation of the cumulative energy 
dissipated by Level Mountain volcanism is presented in 
Figure 9-2. The energy ploted is the sum of three terms 
discussed above. The graph should be read as: 0.979x10*'J 
was dissipated between the onset of volcanism and 6.3 MY 
before present, etc. Two and possibly three distinct stages 
are present. The shield building basaltic volcanism (units 1 
to 4) has the highest rate of energy dissipation, such that 
79% of the total energy expended had occurred by 4.2MY 
before present. The overall rate (which appears to be 
constant) was 1.87x107'J/M.Y. For the plateau surface area 
of 3324 km?, this rate is equivalent to 17.8 mW/m?. The 
closest heat flow measurements in northern British Columbia 
are Dease and Hotailuh which average 73.5 mW/m? (Jessop, 
Souther et al, 1980 in prep.). Presuming the heat flow to be 
Similar at Level Mountain during plateau building times, 
then volcanism represents 24% of the total heat flow. The 


next stage, roughly equal in duration, corresponds to the 
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Figure 9-2. . 
Cumulative energy dissipated by Level Mountain volcanism. 
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bimodal volcanism and stratocone building. Eighty two 
percent of the energy dissipated in this stage is by salic 
lavas. The rate here (also constant), is 0.45x10?'J/MyY, 
which is only 25% of the previous stage. However the surface 
area of the stratocone cap is only 20% of the plateau as a 
whole; expressed in heat flow units the rate for the 
Stratocone stage is 7.15 mW/m?. This is still less than half 
of the shield rate, and only 9.7% of the modern heat flow 
values. 

The final 1.2MY, with its low volume hawaiites and 
alkali basalts, represents another reduction in the rate and 
possibly signifies the extinction of Level Mountain. However 
if the rate of the statocone stage is extrapolated to the 
present, the time period is about right for the release of 
another fifty cubic kilometers of lava. Mountain and for the 
Stikine as a whole. While the preceeding discussion dealt 
with the development of Level Mountain volcanism and the 
total energy dissipated by the process, the following 
section considers the minimum thermal requirements to 
initiate and maintain the volcanism. Basic data for these 
calculations and the results are presented in Table 9-9. 
Assuming the mantle source region to be at the solidus 
temperature, two additional energy requirements must be met: 
(i) the heat of fusion for the melts and (ii) the potential 
energy change associated with the magma transport. The 
energy dissipated through normal heat flow is calculated for 


Level Mountian and, the total volcanic area within the 
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Table 9-9. Energy Requirements of Volcanism 


LAVA Volume (Km?) 

LAVA Mass (10'!?Kg) 
Surface Area (Km?) 
TIME SPAN (my) 

Heat of Fusion (cal/g) 
SOURCE Votume * Km? 
SOURCE MASS 10!?Kg 


Cp U/kg K for 1050°C 
4 phase spinel peridotite 


Modern heat flow mw/m? 
Total heat flux for 7my (J) 
a) H (Jd) 
2) Mgh (J) 
a EVN CASE [ale (Qu) 

Ware 
equivalent heat flow mw/m? 


i aaiein ess ox Gife, (fs) 
available for source heating 


II. Total heat budget for 96° 


heating (uJ) 
equivalent heat flow mw/m’ 


* Source volume has same surface 


Mass conversion assumes mantle 
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Stikine for the last 7MY. For the Stikine this value is 
1.4x10?°J. For Level Mountain this works out to 3% while for 
the total Stikine Volcanic region it is 2.3% of present heat 
EFlow for northern British Columbia. 

Additional energy would be required for pre-heating of 
the mantle source region. For this calculation model source 
region dimensions are estimated from the surface outcrop of 
volcanism and a slab thickness of 4.7km. Source region 
density was assumed to be 3.3g/cc. The heat capacity for a 
four phase spinel periodite was calculated from data in 
Robie, Hemingway and Fisher (1978) to be 1300 J/kg-°K in the 
temperature range of interest. The heat capacity integral 
was solved for a temperature interval by assuming an amount 
of heat potentially available. Less than 40% of a surface 
heat flow value of 73.5 mW/m? probably arises from a steady 
state mantle contribution. A value of 12% of the total heat 
flow is equivalent to 30% of the mantle heat flux. This 
implys that less than 100° of heating can be performed over 
a 7MY period. Obviously heating of this type requires 
geologically long periods of time, or high heat fluxes, or 
small source regions in some combination. For the 12% 
pre-heat budget cited above, the pre-heat requirement is 80% 
and 84% of the total energy requirements for Level Mountain 
and the Stikine respectively. If the energy requirement for 
volcanism is expressed in heat flow units, these heating 
models represent about 15% of the modern heat flow discussed 


above. 
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The low volume alkali magmatism for the Late Cenozoic 
and the Stikine could be initiated and maintained by a 
thermal contribution equivalent to a small percentage change 
in the normal heat flow. If the source region must be heated 
to the peridotite solidus through an interval of a few 
hundred degrees, in the space of a few million years, an 
increase equivalent to 25% of the normal heat flow may be 
required. On the other hand, an increase in available heat 
as low as 2.5% could initiate and maintain alkaline 
volcanism from a source region already at the periodite 
solidus. 

As a final note, the total energy dissipated through 
volcanism in the Stikine can be compared to a change in 
tectonic strain or potential energy. This can be likened to 
a release of mantle strain energy through subsidence of the 
crust over the entire Stikine Volcanic Belt. The area of the 
entire Stikine Volcanic Belt 1S eStimated to be 7.19x10'°m?. 
For 25km of crust of mean density 2800kg/m? the total mass 
of crust under the Stikine is 5.03x10'*kg. Using the 
expression E = mgh to equate the dissipated heat and 
potentialtéenergys(wheres Et] 1030x0227 me=e5e03x108%kg ‘and 
g = 9.81 m/s?) h can be calculated*to be 225m. The efficient 
release of about 0.2km of strain build up over an area the 
Size of the Stikine Volcanic Belt is equivalent to the 
energy released as volcanism. Whether such tectonic strain 
preceeded the Late Cenozoic volcanism or actually was 


responsible is unknown. The release of a geologically small 
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amount of tectonic strain through a melting process could be 
responsible for this type of low volume alkaline volcanism 
without recourse to unusual compositions or events in the 
upper mantle, or large scale lateral geological movements 
such as treated in plate tectonics theory. A model of strain 
release through volcanism is particularly appealing for the 
Intermontane Belt as the strain release through seismicity 


is anomalously low (Milne et al, 1970). 


MODELLING OF A HEAT FLOW AND A GEOTHERMAL GRADIENT FOR LEVEL 
MOUNTAIN AND THE STIKINE 

The intent of this section is threefold, first to 
review existing Cordilleran heat flow data, second to 
calculate a steady state conduction model and third to 
construct a preferred PT model for the geothermal gradient 
under the Stikine and relate it to the origin of the Level 
Mountain volcanics. The steady state conduction model is in 
general, representative of continental geotherms. These 
continental geotherms are characterized by surface heat flow 
in equilibrium with heat flowing into the base of the 
lithosphere, plus a term for internally generated heat from 
radioactive decay (Pollack and Chapman, 1977). 

A review of existing Cordilleran and Canadian heat flow 
measurements is presented in Table 9-10. It is noteworthy 
that the value of 73:.5 mW/m*, used for the Stikine in the 
preceeding section, is virtually identical to the Western 


Canadian average. The Cordilleran values are somewhat higher 
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Table 9-10. A summary of Cordilleran and Other Canadian 
Heat Flow Data 


Gordrllera mW /m? Corrected! Gradient °C/km 
Buckley lake 67 74 
Dease g:1 100 
Hotailuh 56 62 
Penticton* 67 19 34.7 


Canadian, average heat flow by region, uncorrected 


Maritimes (3) Bou5 Eastern Canada 
Quebec (4) 32%2 

Ontario (12) 41.0 Sonte2 
Manitoba (Ae) S305 

Albertaene; (2) 64.1 Western Canada 
NEwWeT? (1) SBea LBYEeE 
Cordillera (4) 7083 


data from (Jessop, Souther et al 1980 in prep.) 
(Jessop and Judge, 1971) 
(Lee and Clark, 1966) 


1. Corrected for the effect of Pleistocene glaciation after 
the measurement of Jessop (1971). 
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than the continental average of 42 mW/m?, but not so high as 
many of the circum-Pacific regions which may have values as 
high as 170 mW/m?. Northern B.C. lies on the northern 
extension of the Western North American thermal high, with a 
regional heat flow of about 70mW/m? (Pollack and Chapman, 
1977). For the Cordilleran thermal anomaly zone, there is a 
sharp eastward gradient of decreasing heat flow (Roy et al, 
1968). Within the thermal anomaly both the crust and mantle 
are hotter than to the east. The corrected shallow 
geothermal gradient deduced for Pentiction (Jessop and 
Judge, 1971) is also high for a continental geotherm. 

A model for heat flow at Level Mountain has been 
constructed assuming a simple three layer geometry. The 
model, in general form, is steady state conductive heat 
transfer through a series of arbitrarily wide parallel 
plates of known thickness and conductivity. The solution to 
Laplace's equation is given by: 

O/At = AT/SUM(1(i)/k(i)). 
The lefthand side of the equation is heat flux per unit area 
per unit time, and the right hand side is the net 
temperature difference across all the layers, divided by the 
sum of the individual layer thicknéss-thermal conductivity 
quotients. Q/At is replaced by a heat flow value and the k 
corresponds to geologically reasonable choices for layer 
compositions. With these values set, the thicknesses and 
temperatures are no longer independent variables so that a 


geothermal gradient can be calculated within acceptable 
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petrological and geophysical constraints. 
The preferred model for steady state heat flow is 81.7 
mW/m? with a two layer crust and one layer upper mantle 


characterized by: 


1 (km) k (W/m°R) 
erust a oO 3.14 
Crust <2 2395 PAM 
upper gmantbee  - \ ---=- ee 


The major limitation of this model is probably the steady 
State assumption. The model and selected solidus/liquidus 
data are presented in figure 9-3. Higher heat flow values 
would cause unrealistically high geothermal gradients that 
do not agree with geophysical or petrogenetic constraints. 
Lower heat flow values (10%) would depress the melting zone 
to depths below 48km which is considered unlikely. 
Variations at this level correspond to lateral temperature 
gradients of less than 0.6°C/km, which could explain the 
narrowness of the volcanic areas in the Intermontane Belt. 
Such small differences be supported by the natural variation 
in thermal conductivity with temperature for peridotites, by 
changes in proportion and composition of phases, or by very 
small amounts of tectonic stress (<1.0km). Caner's (1969) 
estimate for the Curie isotherm depth appears to be too 
great and his 750°C Moho isotherm appears to be too cool. 
The implication of this model is that the Moho under Level 
Mountain is about 25.4km, which would make it one of the 


shallowest values for the Cordillera. Thos model is 
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| Figure 9-3. 
Geothermal gradient for the Stikine for model calculations 
and assumptions as discussed in text. 


A MODEL GEOTHERM FOR THE STIKINE 
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generally consistent with a mantle derivation for the Level 
Mountain lavas formed by the anhydrous partial fusion of 
Spinel peridotite. 

A characteristic of these steady state conduction 
models 1S a steeper geothermal gradient over most of the 
crust than in the upper mantle. For any reasonable choice of 
thermal conductivities, this geothermal gradient behavior 
holds for a wide range of hypothetical heat flow values (73 
to 112 mW/m?), Curie isotherm depths (16 to 18km) and 
crustal thicknesses (25 to 30km). The very shallow 
geothermal gradient is predicted to be lower in the volcanic 
pile as compared to adjacent areas with basement outcrop, or 
immediately beneath the pile in the crust. This is due to 
the high thermal conductivity for lavas. The Intermontane 
Belt in Southern and Central B.C. generally lacks a high 
velocity lid and the LVZ begins at the base of the crust. To 
more accurately fit such a constraint, the gradient in the 
upper mantle would probably be even more gradual. For a 
metasomatized alkali- and silica-rich upper mantle the 
thermal conductivity would be closer to lower crustal values 
and thus support a higher geothermal gradient. Varying the 
model in this direction would also*lend support to a mantle 


derivation of peralkaline melts. 
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Chapter 10. GEOTHERMOMETRY AND GEOBAROMETRY ESTIMATES 


INTRODUCTION 

Spinel lherzolites are a common inclusion type for 
nodule-bearing alkaline lavas of the Intermontane Belt. 
Studies of phase chemistry and trace element concentrations 
have been reported by Littlejohn and Greenwood (1974), 
Fiesinger and Nicholls (1977), Fujii and Scarfe (1981) and 
Fujii et al (1981). Microprobe analyses of coexisting 
minerals from two spinel lherzolites nodules, one from 
Telegraph Creek and one from Grizzly Hill in Central B.C. 
are presented here. The various geothermometers relevant to 
Spinel lherzolites are reviewed here. Equilibration 
temperatures have been estimated for the two Stikine 
localities: Telegraph Creek and Castle Rock (Littlejohn and 
Greenwood, 1974), the Selkirk material (Sinclair et al, 
1977) and Grizzly Hill. The purpose of this is to estimate 
P,T conditions in the upper mantle beneath the Cordillera to 
determine the conditions of alkali basalt genesis for the 
Late Cenozoic, particularly at Level Mountain. 

At equilibrium, the partitioning of Al, Mg, Fe etc. 
between the four phases olivine, orthopyroxene,, 
clinopyroxene and spinel, is a function of bulk composition, 
temperature and pressure. Considerable attention has has 
been given to pyroxenes as geothermometers (Boyd, 1973; 
Wilshue andMJackson, 1975; Presnall, 1975; Stroh, 1976). 


Thermodynamic models may be constructed using tabulated 
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information from mineral thermodynamic studies (Robie et al, 
1978) and experimental data from simplified analogs of the 
natural peridotite system, MgO-Al,0,-Si0,, (MacGregor 1974; 
Fujii and Takahashi, 1978) and the enstatite-diopside join 
(Mori and Green, 1977, 1976; Davis and Boyd, 1966; Lindsley 
and Dixon, -1975) 

Consider the partition of magnesium between coexisting 
orthopyroxene and clinopyroxene 


Mg,Si,0O, = Ma, 51.0% 
Opx Cpx 


At equilibrium, (P,T, composition fixed), 
Opx cpx 
U cow yj 
Mq,Si.,0, Mg, S120, 
the chemical potential of enstatite into both pyroxenes is 
equal (the sum of the u(i) for the reaction is zero). Taking 
the standard state for each component to be equal to the 
pure phase, at some arbitrary temperature, the reaction can 
be rewritten as: 
Opx cpx 
OM BiRDiin ita =SE err RT ola ia 
En En En En 


and the free energy for the reaction is then given by: 


EpxeL Ops 
CAGH).® =5-RieCLn E /a 

Jk En En 
dividing through by RT and recalling the Second Law of 
Thermodynamics: 


AGS = VAN SERS. 


the preceeding expression can be rewritten as 
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to give the basic form of the solid assemblage 
geothermometry expressions. Knowledge of AH° and aS°® for the 
particular reaction in question may come from tabulated 
sources for minerals (such as Robie et al, 1978). 
Alternatively, AH®° and AS° may be estimated from P-T slopes 
of reaction boundaries in experimental phase equilibrium 
Studies using the Clapeyron relation. 
The general form of the solid phase assemblage 

geothermometers is: 

AH°® ASS 

ng eee) am pacar 

RT R 
where K is the equilibrium constant (expression of activity 
products) for the reaction, and AH°® and aS° are the enthalpy 
change (heat of reaction) and entropy change (order at a 
given T and P). As long as the AH°® and AS° are accurately 
known and the simple mixing approximation applies, the 
geothermometers should be concisely expressed as simple 
log-linear expressions ln K = A/T + B, where A and B are 
constants. This is the form generally given in 
geothermometry papers (Wells, 1977; Mori, 1977; Herzberg and 
Chapman, 1976) but it is important to keep in mind that A 
and B are not just arbitrarily selected numbers chosen to 


best fit a particular collection of data, but that they have 
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fundamental thermodynamic significance to them. Where the 
ideal mixing assumptions break down, or data on simple 
experimental fail to fit the more complex natural 
assemblages, empirical formulations have been made using 
polynominal expressions and adding additional terms, such as 
Mori (1977) or Wells (1977). It should be remembered that 
these extra terms are really without physical meaning and 
have only been added to force a mathematical fit toa 
collection of data. The need for such modifying terms only 
underlines the inaccuracies in our knowledge of AH®°, AS° and 
in the assumed solution models. Herzberg and Chapman (1976) 
point out their fits imply aS >8 cal/mol° which does not 
seem high in this temperature range. Mori's values (1977) 
imply even higher a S (> 9.4 cal/mol®). This probably 
demonstrates the oversimplification of the mixing model, 
rather than inaccuracy in the thermodynamic or experimental 
data. The familiar thermodynamic form of the geothermometer 
can be rewritten aS an equivalent polynomial in In K so that 
all geothermometers can be expressed in a common form: 
10270 Ke=e avin kk) a2 pune) tc 
where c 10CKS AH - 


b 104R/ AR 
and a is an additional term with*no thermodynamic meaning. 


The state of the art in geothermometry of spinel lherzolites 
is summarized in Table 10-1. In addition to the magnesium 
partition between pyroxenes used in the derivation above, 


(reaction 1, Table 10-1) the two most accurate 
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Table 10-1. GEOTHERMOMETRY OF SPINEL LHERZOLITES 


ASSEMBLAGE : Forsteritic Olivine - ol 
Crzalepiopside — cpx 

Magnesian Orthopyroxene - opx 

(Mg,Fe,Ai,Cr) Spinel phase FSH 


REACTIONS: AT EQUILIBRIUM TO DEFINE P,T 


CDK Lops 
1)=Mg>551,0; = Mg3S12.0¢¢6 (k1) elifa Ja ) 
Opx Cpx Prt En En 
2) CaMgSi.,0% # MgAl,O% = CadAlgS1OpetwMgsgS20.5 
opx Sp cpx ol 
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Table’ 10-1... Ccontinued) 


ACTIVITY DEFINITIONS' 


PEAGOxenes MMC Ala Coy Vow bu jhe), ber" Mg** »Ni?')), 
M2(¢Mm4 39Ca 7", Na Diba afFe? t,Mg?* ) 


Cpx M2 M1 In all cases Al** is partitioned 
a = X werx between the tetrahedral site and M1 
En Mg Mg octahedral site on the basis of 
Silica content in the pyroxene 
cpx M2 M1 analysis assuming the tetrahedral 
a =X .X Site is completely filled by 
Di Ca Mg Si** and Al?*. After assigning 
the site specific cations Mg’*, 
cpx M2 M1 Fe?* are partitioned according to 
a = x Se their proportions between M1 and M2, 
CaTs Ca Al assuming ideal mixing. 
Opx M2 M1 
a = & see 
En Mg Mg 
Opx M2 M1 
a = x Sars 
Mats Mg Al 


(Olivine) 


ol ol 
a = (xX) i) 7eerdeals solution iM) .S10% 

Fo Fo 

(Spinel) 

Sp A B 
a = ox (X ea)etiowmiiormuvaeABeOe normal spinel 
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Calibrated Geothermometers: 
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Linear or polynomial regression fits of the form: 
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Activity definitions conform to Nicholls and Carmichael 
(1972), Bacon and Carmichael (1973) and reference cited 


below 


Herzberg and Chapman (1976), data from experimental and 
natural spinel peridotites 


Mori (1977), experimental data base 
Mori (1977), after MacGregor (1974) 


Mori (1977), after Obata (1976) 
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geothermometers are based on magnesium and aluminum 
partitioning between coexisting clinopyroxene spinel and 
olivine (reaction 2, Table 10-1) or aluminum and magnesium 
partitioning between coexisting orthopyroxene, spinel and 
olivine! (reaction 3, Table 10-1). Reactions 2 and 3, 
involving the solubilities of alumina in pyroxenes, are not 
yet tested for Fe- and Cr-bearing systems. Application of 
theses geothermometers to the Cr bearing system may give too- 
high a temperature estimate (Fujii and Takahashi, 1978). The 
definitions of activities are consistent with the 
geothermometry papers previously cited. Activity expressions 
are calculated for ideal structural formulae and using 
concepts of equivalent sites. For example the expression 

ol M1 M2 M 

Bea = Ree. nie hee ON for M1 = M2 

Fo Mg Mg Mg 
Shows the activity of the forsterite component in olivine 
assuming a random mixing for iron and magnesium between two 
octahedral sites. There is ambiguity in calculating 
Structural formulae, particularly for pyroxenes, and several 
methods have been proposed (Cawthorn and Collerston, 1974; 
Fleet, 1974a,b). They have shown for pyroxenes of 
intermediate iron content, that random partitioning of Fe’* 
between the two octahedral sites iS not a very accurate 
assumption, because Fe?* partitions in favour of Mg’* into 
the M2 site. Magnesium and aluminum rarely exist in spinel 
analyses in stoichiometric proportions unbalanced by other 
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laboratory studies. Additional difficulties may be 
encountered when fitting natural spinel lherzolites which 
have many additional components such as Ni, Mg, Na, Cr, Ti, 
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REVIEW OF PREVIOUS DATA FOR THE INTERMONTANE BELT 

Fiesinger (1975) and Fiesinger and Nicholls (1977) 
performed geothermometry calculations both on solid 
assemblages and silica activities with coexisting lavas for 
spinel lherzolite and gabbroic xenoliths from Late Cenozoic 
lavas of the Wells Grey - Quesnel Highlands area B.C. 
Calculated equilibration temperatures ranged from 1050° to 
1650°C with a total range of equilibration pressures from13 
to 50 kbar. The larger values for P and T here lie outside 
the spinel lherzolite stability field even though no garnet 
lherzolites are reported. There is too much variation in 
calculated P and T to be explained by any fractionation 
processes at either high or low pressure or to be explained 
by the inaccuracy in the assumptions and calculations of the 
geothermometers. The implication is that the nodules have 
been sampled from different depths and that the geotherm 
varies both in time and space for the Intermontane Belt. 
Similar conclusions have been obtained by Fujii et al (1981) 
Using data tlorietwo Glocalities.sirom (central sand hsouthernrB .G% 
Takomkane (Boss Mountain) and Clark Creek (West Kettle 
River). Littlejohn and Greenwood's (1974) data for Jacques 
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Stikine Volcanic Belt) indicate high and variable 
equilibration temperatures according to the olivine-spinel 
geothermometer of Jackson (1969). These estimates should be 
considered with some doubt however, due to possible 
subsolidus reequilibration of olivine and spinel and the 
assumption of ideal solution in the spinels. Sinclair et al 
(1977) also calculated high equilibration temperatures for 
nodules from the Selkirk cone. 

Chemical analyses of four phases were made by electron 
microprobe on polished block mounts of nodules from Grizzly 
Hill and Telegraph Creek. Data and structural formulae are 
reported in table 10-2. All of these analyses gave high 
totals between 101 and 106% due to problems with an oil film 
on the beryllium window on the microprobe at the time the 
Gata were collected. Despite the high totals, the recast 
analyses compare favorably with other mantle inclusions 
reported from the literature including samples from the 
Cordillera (Littlejohn and Greenwood, 1974). 

Activities have been calculated, from the structural 
formula and site occupancy restrictions as summarized in 
Table 10-1, for these nodules and for the four Castle Rock 
samples from Littlejohn and Greenwood (1974), and the two 
samples from the Selkirk Cone (Sinclair et al, 1978). 
Temperatures have been calculated for the geothermometers of 
Herzberg and Chapman (1976), Mori (1977) and Wells (1977) 
for all four localities, see table 10-3. The Castle Rock and 


Selkirk temperatures are derived from averages based on all 
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Table 10-2. CHEMICAL ANALYSES OF FOUR PHASES FROM INCLUSIONS 


A. OLIVINES 


TCL GHL 
1004-41 
Si0, 40.58 40.48 TCL 1004-41 also contained 
interstitial plagioclase <1% 
FeO 10.95 is25 of An(55)Ab(43)0Or(2) and mono 
sulfide solid solution phase 
MnO 09 oS approximately Py(64)Pn(36) 
having detectable Ni,Co,Zn and 
NiO 40 247 micro-inclusions of sphere 
GHL additionally contained <1% of 
MgO 47.90 47.67 interstitial alkali feldspar 
of An(0.7)Ab(41.0)0r(58.3) and 
CaO 09 0.0 a Similar mono-sulfide solid 
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STRUCTURAL FORMULAE PER 4 OXYGENS 


SP 6). 1. 00#81 20 7000'4.0 
Fe?* .226 ES 
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Table 10-2. 


B. Pyroxenes 


Cpx 
2G GHL 
1004-44 
S203 53656 53.49 
P10, Ce28 O35 
Ad 0; 5264 Seo 
ae AO O25 0267 
MeO; OF 6), Oeh1 
FeO 2379 Ze to 
MnO 04.0 0.0 
MgO 14.25 14.28 
CaO 22223 22.40 
NiO 0.0 0.04 
Na,0O 0.69 0.63 
BaO 0.06 0.0 
99.99 100200 
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Cr>: 021 
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Fe?* .084|1.946 
Mn? * —Ga 
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CO 001 
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Table 10-2. (continued) 


C. Spinels 

TEL GHL 

1004-41 

04 0.09 O20 
Ale Ose 54.52 59.0 
Cree Hioo 8.99 
Ver 0.0 G..0 
Fe,0; 4,46 8.99 
FeO Teta | 9.28 
MnO OZ 0.09 
MgO 2s Se 20.62 
NiO 0.43 0.45 
zno Gero 0.05 
CaO 0.14 ae 


100.00 100.02 


STRUCTURAL FORMULA PER 32 OXYGENS 
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Table 10-3. Calculated Equilibration Temperatures (T°C) for 


Spinel Lherzolites from the Intermontane Beit. 


(1) (2) (3) (4) 


SELKIRK TELEGRAPH CASTLE GRIZZLY 
CONE , YUKON CREEK EB. ©. ROCK? Be Ci iwHILetWBeC. 
72 Kbar 
*K , 12 0 10393 2 Fi 10414 
*K » 1304 1366 1261 1S 
76 Kbar 
*K, 1151 1066 1308 1069 
+Ku 1126 1010 1304 HOLS 
*K, 1323 1381 1264 1336 
aa sts 1246 12993 1208 1258 
STK, Pega 1246 1324 1146 
#K 5 13093 1267 1388 1alse7. 
O-40 Kbar 
#K , 962 878 1095 878 


(1) 


(2) 


(3) 


Sinclair, Tempelman-Kluit, Medaris (1977) for mineral 
analyses of 2 Lherzolite nodules in Pleistocene cinder 
cone, central Yukon 

Stikine region: Lherzolite nodule from Miocene or 
younger basalt flow at Telegraph Creek, B.C. Sample TCL 
1004-41 from H. Greenwood 

Littlejohn and Greenwood (1974) for mineral analyses 
from 4 Lherzolite Nodules in Quaternary AOB breccia at 
Castle Rock (Stikine Region) on northern edge of 
Klastline Plateau 

Central B.C., Lherzolite nodule from Miocene basalt 
dyke at Grizzly Hill Hamilton (1979) 

Wells (1977) 

LWOF ATO Kee 2153626 lint kite 52.008 6.1324. Fe/Pe'+\ Mg: ops) 
Of the natural compositions he tests this on, 
lherzolites and peridotites give the largest T°K 
discrepancies compared to other mineral geothermometers. 
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of the reported mineral data. The calculated temperatures 
fall in a tighter range than for previously reported values. 
Taking the Telegraph Creek, Castle Rock and Selkirk values 
as representative of the Upper Mantle beneath the Northern 
Cordillera in the vicinity of Level Mountain, the three 
mineral thermometers (k.,k;) give values in the range 1200° 
to 1400°C at 12 to 16 kbar, while the pyroxene thermometers 
give lower estimates of 875° to 1300°C. Using a surface 
temperature estimate (5°C), the average of all 12 kbar 
estimates (1219°C), and the average of all 16 kbar estimates 
(1228°C) for each locality, results in geothermal gradients 
in the range of 25°/km (82.4°/kbar) down to the intersection 
with the lherzolite solidus, see Figure 10-1. Owing to the 
Cr content of these natural samples, this geothermal 
gradient estimate is probably a maximum. For the simple four 
component experimental analog of the lherzolite system 
(CMAS)'this intersection would be at 15.8 kbar (51.8 km). 
Considering the effect of other components such as alkalis, 
iron, titanium, phosphorous and volatiles in depressing the 
solidus, this intersection is likely to be even shallower 
(see pyrolite solidus of Green and Ringwood, 1967). If some 
of the higher geothermometry estimates are accurate, the 
geotherm could be as steep as 31°/km with a solidus 
intersection as shallow as 41km (12.5 kbar). On the basis of 
trace element data, all of these spinel lherzolites appear 
to be high temperature refractory residua from previous 
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Figure 10-1. 
P-T plot showing phase fields and solidus for mantle 
peridotite along with liquidus data for basaltic lavas and 
range of PT estimates for lherzolites from northern 
Cordillera. 
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partial melting events and they are too depleted to be 
source rocks for normal basalt genesis. 

For comparison with natural rock systems the melting 
relations (solidus and liquidus) of a basanite are shown in 
Figure 10-1. The baSanite comes from Mount Shadwell, 
Australia cArculus, 1975).. The intersection of the basanite 
solidus with the hypothesized Stikine geotherm lies in the 
Stapmlit yervelaror obi aCow teSpict Plo + Ne andwethe i] rquidus i 
geotherm intersection is in the field of ol + L, according 
to the phase relations of Arculus (1975). 

From the work of DePaolo (1979) the intersection of the 
liquidus of an AOB having log a S$i0, = -0.4 with the 
anhydrous pyrolite solidus, is at 13.5kbar (44.5km). This 
point marking the alkali olivine basalt depth of origin lies 
within the field of the northern Cordilleran geothermometry 
estimates and slightly to the high temperature side of the 
hypothesized Stikine geotherm. Using these data the melting 
region for generating the Stikine basic lavas can be 


approximated to lie between 36 and 50 km depth. 


GEOBAROMETRY AND GEOTHERMOMETRY FROM THE ACTIVITIES OF 
SILICA AND ALUMINA IN LEVEL MOUNTAIN LAVAS 

Expressions for the variation of log aSiO, into the 
melt with temperature have been published by Carmichael et 
al (1970) and for temperature and pressure by Nicholls et al 
(1971), Nicholls and Carmichael (1972) and Bacon and 


Carmichael (1973). The basic assumption is that for the 
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coexistence of minerals and silicate melts at equilibrium, 
the compositions of the minerals define a buffer on the 


activity of SiO, into the melt. Thermodynamically this can 


be stated: 
melt melt minerals minerals 
ue + RT -ln’a ew be FeRT Anca 
Si Or SiOs S10, SiO; 


and for more than one solid phase the right hand size would 
contain additional activity terms. If the chemical potential 
of silica in the standard state is taken to be that of SiO, 
glass for both melt and crystals, then the activity terms 
are equal. The activity of silica into the melt (or solids) 
may be defined by: 


aly it ye ple $ ------- (Phe + logrk 


2230 SRT, 2%, 30 3RT 


melt AG® AV? 
[10g a 
pe? 


Sigs 
# where ACp = 0 and K is the equilibrium constant for the 
reaction in question, defined by stoichiometry, and the 
activity terms are for pure mineral components into actual 
crystal phases, after Bacon and Carmichael (1973). Examples 
of such reactions that buffer silica are: 
Si0y (melt) i=4S 10308 tOuartz) 

andl Mg¢.$3a 03) @, SiO5\( men tavi=H2Mg64 Os 
after Nicholls et al (1971). Note that in using real mineral 
compositions, the activity depends on both the solution 
model assumed and the stoichiometry. The thermodynamic form 
of the silica activity stated above gives rise to a unique 


expression for each reaction of the general form: 
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melt 
[109 a | eT eB er Cmr en Pen)i) at log ak 
Si Gael Rar 


In keeping the same form, it can be seen that A and B come 
from the free energy term, while C comes from the volume 
change. Values of A,B and C have been calculated from 
existing thermodynamic data for the various silica buffering 
reactions. In the absence of a mineral-melt reaction which 
defines (buffers) the activity of SiO, into the melt, the 


following expression is used 


melt melt W(Si02) Va 
[10g a | = [10g Xx | z 
Det 2. 30R8 


S$i0, SiO, ae 
after Bacon and Carmichael (1973). The V term is evaluated 
uSing the partial molar volume data of Bottinga and Weill 
(1970). The other term is a function of composition 
melt 
where X is the mole fraction of silica into the melt, 
SiG. 
and the composition function is defined: 
melt melt 
$(S10,)/T = log E ait 
S207 S10, 
The function $(Si0O,) may be evaluated at any single 
temperature if the activity of Si0O,° into the melt is known 
Nicholls et al (1971) and Bacon and Carmichael (1973) have 
evaluated the activity of $10, at the quench temperature of 
a lava. For this purpose they used T quench and log f0O, from 


Buddington and Lindsley's (1964) Fe-Ti oxide data, and 


equilibrium groundmass mineral compositions for olivine and 
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titanomagnetite with the fayalite-magnetite-silica buffer. 
The buffer expression is: 

melt ol TiMt 
rogca = A/Tat Bet og a Pt SLOG fO, 2) 2/3. Log ra 

$10, Fa Mt 
Note the reciprocal dependance on T and log f£0,. A small 
error in estimating T and/or log f0,, while it changes the 
one bar activity into the melt, does not have as great an 
effect on estimating $(Si0O,) for the lava. The confidence of 
X(SiO,) depends only on the accuracy of the bulk chemical 
analysis. 

Alumina activities are evaluated in a parallel fashion 
for mineral reactions that buffer Al.,0,. For the unbuffered 
case, log aAl,0O,; into the melt must be evaluated at some 
known temperature to get at MAl1,0;. In basic and 


intermediate lavas a useful defining reaction is: 


CaMoS107 (cpw) #2172551 05 (mel ayer eAlO, (melt) 


= CaAl,Si07,¢plag) + 1/2: Mg.Si0,(01) 


where the equilibrium groundmass assemblage is used and log 


aSiO, into the melt must be known. The expressions 


melt pl ery 2 cpx melt 
log a = A/T+Btlog E (ae ua lr i/2log a 
PelnsOs An Fo Di GEO, 
melt melt 
and QA1,0; = [109 a Skog a2 |r 
A lesOrs AOR 


are evaluated using published values for A,B (Bacon and 
Carmichael, 1973) and compositions data for a specific 


basalt and its constituent mineral phases. 
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The variation of mineral-melt buffers for aSiO, and 
aAl.O; are shown with temperature and pressure in figure 
10-2. These diagrams have been constructed using the data of 
Bacon and Carmichael (1973), Nicholls and Carmichael (1972) 
and Nicholls et al (1971). Judging from common phenocryst 
assemblages, these curves likely span the aSiO, and aAl,O, 
ranges commonly encountered for igneous melts from basalts 
to rhyolites. The unbuffered curves shown for basalt are for 
examples from Level Mountain. While the convergence of the 
common silica buffers at high temperature limits the 
applications of log aSiO, into the melt as a geothermometer, 
the strong pressure dependence for the mineral-melt buffers 
compared to unbuffered lavas lends itself to use as a 
geobarometer. Log aSiO, varies inversely with pressure, 
while log aAl.,0; varies directly. The use of silica and 
alumina buffers together improves the precision of pressure 
estimates. 

Having defined the unbuffered variation in log aSi0, 
and log aAl,O, intelithe meltwas a function "of ¥Ptand T for a 
given lava, the expressions for both log aSi0O, and log 
aAl,O, may be solved simultaneously with reactions that 
define the activities of those components in the melt for 
mineral buffers, using assumed or analysed phenocryst or 
mantle assemblages, to give unique curves in PT space. As an 
example, Bacon and Carmichael (1973) solve the 
forsterite-enstatite buffer for SiO, (defined previously) 


and the alumina buffer given by 
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Figure 10-2, 
Calculated alumina and silica activity mineral-melt buffer 
Curves. Pure phases unless otherwise shown. Data from 
Nicholls and Carmichael (1972) and Carmichael (1973). 
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Mg2Si0s(ol) a -Al70,(melt) s=eMgAl,O4(¢sp) + MgSi0, (opx) 
with the unbuffered log activity expressions for the lava to 
give an assumed P and T of equilibration with mantle 
lherzolite. From the joint solution of the expressions for 
any mineral assemblage buffer, the general P-T curve is 


given by: 


where T is in °K, P in bars, A,B and K refer to the 
particular mineral-melt buffer and D+tET comes from the 
difference in the C coefficients for unbuffered melt and the 
chosen buffer reaction. In the case of the 
forsterite-enstatite buffer this expression is: 

(0-597 -—elogE( SiOn mr Log ik)it) = (OS1O7 8431034) 

ORO See hi tn 47 ex 105" 1 

where log K depends on the composition of olivine and 
pyroxene in the lherzolite. Similarily the 
olivine-orthopyroxene-Sspinel alumina buffer and the log 
alumina activity expression for the unbuffered melt can be 
solved to give: 


(2Gsi2 tee loge Kin alOG mn Ure Og) te (OAS OcmerarG loki) 


where log K is given by: 

logike= logsatSpmimeisyo) +thlog a(n iimopx) + log a\ho anol) 
The equilibrium P-T curve for the lherzolite can also be 
‘delined @by#thesceact son. 
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and ~hogikr="=2156/T + 07069 "0.00778 (P-1) 

The geobarometry of salic lavas presented in Nicholls 
et al (1971) depends on the simultaneous solution of two 
mineral-melt silica buffers. The fayalite-magnetite and 
quartz-Silica buffers were used for a variety of acid lavas 
from different localities including a pantellerite from 
Pantelleria with _P- eq. -> 7.3 kbar. 

Following the method described above, similar 
geobarometry calculations have been attempted for a few 
representative lavas from Level Mountain. The data used in 
these calculations and the results are presented in table 
10-4. The quench temperatures inferred here are generally 
less than or equal to solidus estimates for Level Mountain 
lavas of similar composition. The oxygen fugacity inferred 
is slightly greater than or equal to that of the melting 
experiments. buffers for FMQ and MW by gas mixing. 

While no lherzolite inclusions were found in the Level 
Mountain lavas proper, lherzolites have been found in 
Similar lavasstrom-the Selkack}Cone (Sinclain etsaks 1978), 
Castle Rock (Littlejohn and Greenwood, 1974) and at 
Telegraph Creek. The mineral compositions from those 
localities have been used with Bacon and Carmichael's data 
£or 

CaAl,SHO, * Mg;Si0, = €aMgSi,02 * MgAl.0, 
to give the following lines in PT space for the spinel 
lherzolites from the Northern Cordillera (for T in °K and P 
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Telegraph Greek: «sPo= 275.494 (B= 272122 
Castle Rock: Be= oT Goal ehesas/ 2122 
Selkirk: PasiG reste. eAZse 7122 


The Castle Rock assemblage gives a very Similar slope 
(AP/AT) to the San Quentin» lherzolites.of ‘Bacon and 
Carmichael (1973), while the slopes for Telegraph Creek and 
Selkirk are slightly less. Graphical presentation of the P-T 
solutions for a basalt from the base of the Level Mountain 
plateau (unit 1) and a hawaiite from the central stratocone 
(unit 8) are given in Figure 10-3. Not only do these 
examples span most of the time of the Level Mountain 
volcanism, but they also span the range of depth estimates 
for the basaltic lavas. Two of the hawaiites presented in 
the table do not have real solutions (lying in positive P-T 
Space) for lava-lherzolites or lava-phenocrysts. Those two 
hawaiites are assumed to be derived by substantial 
fractional crystallization from more primitive compositions, 
and in their final erupted compositions are likely to be out 
of equilibrium, with respect to silica and alumina activity 
for their phenocrysts and their presumed mantle source. As 
with the San Quintin lavas presented by Bacon and 
Carmichael, there may be some problems with the calculated 
pressures and temperatures. For the Level Mountain basaltic 
lavas all of the pressures fall into the spinel lherzolite 
Stability field, but the temperatures seem too low for 
genesis of basaltic magmas. The pressures are similar to 


those calculated for other alkaline lavas from the. 
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Figure 10-3. 
Calculations Ofer -T equilibration conditions for spinel 
lherzolites from the northern Cordillera and two Level 
Mountain lavas. Upper set of 3 subparallel curves are for 
spinel lherzolites alone, assuming the reaction 
CaTs+Fo=Dit+Sp. The plot on the left for PAX, from the lower 
plateau, shows the range of equilibration pressures for the 
unbuffered melt and lherzolite minerals, using all 3 
lherzolites. Points are the intersection of Fo+Si0O,=En with 
FotAl,0O,=SptEn. An estimate of pressure of equilibration for 
unbuffered melt and phenocrysts 1S given by the intersection 
of Jd+Si0,=Ab for the plagioclase with Di+Si0,+A1,0,=An+Fo 
for the pyroxene. The plot on the right is for stratocone 
hawaiite 8/25-50/6397. The pressure estimate for melt 
equilibration with lherzolite minerals is lower. The lower 
of the pairs of curves for phenocryst buffers (with negative 
Slope) is for CaTs+Si0O,=An other reactions as before. T= 
Telegraph Creek, S = Selkirk, CR = Castle Rock. 
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Intermontane Belt (Fiesinger and Nicholls, 1977); Nicholls 
et al, 1981 in prep.). The pressures calculated for three 
Salic lavas from Level Mountain are all in the range from 
6.0 to 6.5 kbar. This corresponds to the average depths 
calculated for diverse salic lavas in Nicholls et al (1971). 
The uncertainty associated with the depth (pressure) 
estimates for Level Mountain salic lavas is far greater than 
for the basalts. The two reasons for this are that the 
fayalites for these lavas were in reaction 
(petrographically) with their respective groundmass 
compositions (as judged by riebeckite or sodic pyroxene 
overgrowths), and the mineral buffers used to calculate 
depths were only assumed to be in equilibrium with said 
lavas, and this may not have been the case. The salient 
features of these calculations for Level Mountain are 
1. that most of the basalts could have been in equilibrium 
with a spinel lherzolite mantle at pressures in excess 
of<9.4 kbar 
2. that phenocryst assemblages from most basalts imply 
crystallization in the 4 to 6 kbar region, which could 
indicate the existence of crustal level magma chambers 
and 
3. that salic lavas could have come from the lower crust or 
upper mantle according to the silica activity 
geobarometry calculations for three trachytes. 
A summary of the geobarometry calculations for Level 


Mountain lavas is presented in Figure 10-4, with pressure 
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presented versus elevation (eruptive sequence) of the flows. 
The implication is that initial plateau volcanism originated 
at pressures from 18.5 to 21.0 kbar and subsequent lavas 
were derived from shallower depths. There is an indication 
of an upward migrating melting zone with time. This would be 
consistent with expectation from experimental studies on the 
Spinel lherzolite system (Takahashi and Kushiro, 1981; 
Presnall, 1980; Ito and Kennedy, 1967). Following the 
reasoning of Presnall (1980), for a CO, bearing upper 
mantle, a rising geothermal gradient would intersect the 
lherzolite solidus cusp in the region of 20 kbar and remain 
there until one phase was exhausted (generally 
clinopyroxene). Subsequently the geotherm is constrained by 
the geometry of the spinel lherzolite solidus and a period 
of rising depth of ‘meltingucould occur in a short’ time span, 
until the geotherm intersects the 10 kbar cusp on the 
lherzolite solidus. This sequence of events could have 
occurred at Level Mountain with the second cusp being 
reached by the time of extrusion for the basalts of the 
uppermost plateau (unit 4). The lavas of the stratocone 
could represent the further upward migration of the geotherm 
and melting zone into the base of the crust, or by 
generation of primary magmas at the 10 kbar cusp and further 
evolution in crustal level magma chambers. 

Takahashi and Kushiro (1981) have found a partial melt 
of "phonolite" composition to exist in equilibrium with Fe- 


and K-rich spinel lherzolite (Salt Lake Crater, Oahu) at 
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Figure 10-4, 
Geobarometry and geothermometry estimates for equilibration 
of Level Mountain basalts with spinel lherzolite using 
(Fo-En-Si0,) and (Fo-En-Sp-Al,o,), imply a rising geothermal 
gradient and a rising zone of partial melting with time. 
Phenocryst equilibration conditions range from the shallow 
upper mantle to the lower crust. Trachyte estimates from 
(Fa-Mt-SiO.) and (Q-SiO,) imply equilibration in the lower 
crust. 
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temperatures as much as 200° below the apparent lherzolite 
solidus (anhydrous). These melts contain in excess of 9.0% 
by weight K.,O and are effectively peralkaline. This is 
reminiscient of the bimodal melting hypothesis of Yoder 
(1976). As the stratocone lavas from Level Mountain include 
phonolites and are dominantly peralkaline, the possibility 
still remains that they could be directly generated at 
mantle depths. Takahashi and Kushiro (1981) also performed a 
variety of "pyrolite" type experiments involving basaltic 
glass sandwiched between peridotites (spinel lherzolites). 
They found that plagioclase - pyroxene gabbros crystallized 
below the lherzolite solidus and were apparently stable in 
the PT field for spinel lherzolite. Should significant 
portions of such alkaline gabbros remain behind in the upper 
mantle they could contribute to the mantle differentiation 
and to compositional variation of subsequent partial melting 
episodes . The existence of gabbroic compositions at mantle 
depths could possibly account for the peralkaline salic 


lavas cf the stratocone stage. 


fa oat. antaaete wt + ib hdl 
a) set OK A Ade: at ann au Rta dl 
ee ae wad tr stipe ve lbp 

yp le e Me a On e bo 4 4 
erat + Lent olga ods 0 ins 98 
Site, Pay aioe owed anidarararte pik. ae te 


= 
ted? 


YS) RAD Ls vit rennet #58 bn aedifo 

hie et bene yiey lie > gts ‘yeas, ee # 

‘we 6 LE os baie Bats baat fy jie " stp Mss 

PRET hare | me tart KG: ous ‘ spe a a aes 

‘ities sin 3 idee ngonged toa 

40 fis Sve ki RA a r tpahontlgtag hats oy or 

dant ris Sees pas ean brew nd Sesto i 

aa Lane, iad Perc. hdntige set legit i ne 

ego 6 an» vsuelle ni eee or ae st: ‘tt rf wiht anus oa 
| 47 HEL, flat ‘teat me i 4 seeks lelnioe wi a. 

Pa hect | is tippwedie: Bis ne eatingy hag been . 


SL 706M 25 apes 2 RORNED arouiawe! Ye aphsas pt ae 
oni Led fara aut 254 aie sthetiong stu 


re 


CHAPTER 11. A PALEOMAGNETIC STUDY OF THE LEVEL MOUNTAIN 


VOLCANICS 


INTRODUCTION 

Paleomagnetic studies have been conducted both on Mount 
Edziza (Souther and Symons, 1974) and on Level Mountain 
(Evans and Hamilton, 1980). At Level Mountain fifty cooling 
units were collected from two stratigraphic Bec e one Ma auee 
11-1. The composite thickness of the proximal and distal 
sections 1S greater than 1000 meters. Fourteen magnetic 
polarity zones are found to be present as well as several 
low latitude virtual geomagnetic poles (VGP's) which 
probably represent transitions or excursions. A combination 
of magnetic data, stratigraphy and fission track dates 
indicates that the time span of Level Mountain volcanism 
extends from Epoch VII (greater than 6MY) to the Brunhes 
Epoch (less than 1MY). This will be shown to be similar to 
the results from Edziza. 

Most sites carry a stable remanence (TRM), which can be 
determined with high precision. Fisher's precision parameter 
k generally exceeds 100. The two stratigraphic sections 
include several low latitude VGP's that probably represent 
polarity transitions or geomagnetic excursions. Filtering 
out the sites with low-latitude VGP's (N=5) and high 
internal scatter (N=7) leaves thirty-eight site means, 
twenty-four of which have normal polarity and fourteen of 


which are reversed. The average normal direction does not 
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Figure dit=1, 
Location of the two paleomagnetic sections are shown on the 
geological sketch map of the Stikine Plateau: showing Level 
Mountain range and location of related Late Cenozoic 
volcanics. 
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differ significantly from the inverted reverse direction, 
and the overall mean yields a paleopole at 85°N lat. 256°E 
long., with an A,s of 7°. The actual mean paleopole is 
slightly near-sided in terms of the common site longitude 
convention, which fails to support Wilson's eccentric dipole 
model (Wilson, 1972). Instead near-sidedness seems to be a 
persistent characteristic of paleopoles obtained from 
Western North America, suggesting the presence of a long 
term non-zonal anomaly. 

The paleopole position for the best thirty-eight sites 
at Level Mountain includes the present spin axis as does the 
Edziza paleopole and the paleopole for the combined Stikine 
data sets. Since the circles of confidence for these 
paleopoles include the spin axis, any Late Cenozoic tectonic 
movements for the Stikine Terrane would appear to be small. 
The primary purpose of this study was to establish a 
magnetic stratigraphy for Level Mountain in hopes of 
correlating to the previously published magnetic 
Stratigraphy of Edziza (Souther and Symons, 1974). An 
additional incentive was to seek information concerning the 
overall morphology of the geomagnetic field and its change 
with time. In this regard, the Late Cenozoic age and high 
latitude location of Level Mountain are particularly 
important. The rock magnetic properties of the samples were 
also investigated to document possible correlations between 


magnetic and petrological characteristics. 
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METHODOLOGY FOR THE PALEOMAGNETIC STUDY OF LEVEL MOUNTAIN 
Field Procedures. Two stratigraphic sections of twenty-five 
flows each were collected. Five to seven cores were sampled, 
Prom each Of the “S0 -flow-cooling “units “for a =totalyof- 252 
cores. With three exceptions all of the flows were collected 
from contiguous stratigraphic sections. In these cases PBT, 
PBX, and PAF either exposure or access required that flows 
be projected "line of sight" back into the main section. 

The PA section (distal and stratigraphically lower) was 
collected from the east facing canyon wall of the Little 
Tahltan River on the Southern Plateau margin at 
approximately 58.25°N 131.50°W, Plate 11-1. The section 
crosses 280 meters of vertical relief which spans the first 
four eruptive units (1-4) encompassing the plateau building 
Stage at Level Mountain. 

The PB section (proximal and stratigraphically higher) 
was collected from the South face of Meszah Peak and the 
south trending hanging valley which drains into the 
Kakuchuya Creek at approximately 58.49°N and 131.48°W, see 
plate 11-2. This section crosses 700m of vertical relief 
Spanning the stratocone eruptive units 5 to 9. There is a 
134m elevation gap between the two stratigraphic sections, 
which are separated by a lateral distance of 23km. According 
to field mapping there is no stratigraphic overlap between 
the two sections, but little or no section is inferred to be 
missing. The preferred interpretation is that the 134km 


elevation gap reflects the slope of the plateau surface 
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Plate wis. 
The PA section through the basaltic flows of the Level 
Mountain was collected from the canyon wall of the Little 
Tahltan river, site PAA to upper right, site PAY lowest dark 
flow in stream notch above tree line. Map unit 1 extends 
from base of section to top of dark flows low on cliff. Map 
unit 2 is lighter coloured with thin flows. Map unit 3 
comprises dark flows and tuffs to top of notch. Map unit 4 
rests on a thick light coloured tuff unit that lies along 
the top, ofthe clip, 
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Platemric2s 
View to SW along Hanging Valley from summit of Meszah Peak. 
The PB section through the stratocone followed the creek, 


down the hanging valley into the lower valley of Kakuchuya 
Creek. 
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BO Sas? )t 

Laboratory Methods. The general measurement techniques to 
obtain the remanence vector for rock samples have been 
reviewed by McElhinny (1973) and Cox et al (1967). Pilot 
measurements of the magnetic vectors were made on a 
Princeton Applied Research Spinner magnetometer model SM-1 
interfaced to a Wang 600 programmable calculator. The 
Stability of one specimen from each flow was tested by 
routine alternating field (AF) techniques, mostly in peak 
freldsvof 265r5)wi0s 20.4404vand80 mle iThetresults were 
used to identify regions of directional stability, and thus 
to select appropriate treatments for the remaining 
specimens. The bulk of a lab work was conducted using the 
Schonstedt SSM-1 interfaced to the Wang. The entire study 
involved about 1100 remanence measurements (presented in 
Appendix 3). Only one specimen from each core wasS measured 
as it has been demonstrated by Doell and Cox (1963, 1965) 
that very little additional accuracy is gained by multiple 
Sampling. Each specimen was measured for natural remanent 
magnetization (NRM) and stepwise demagnetized to find a best 
estimate of D and I for each site. The stepwise 
demagnetization was performed in an alternating field (AF) 
apparatus designed and built by Murthy (1969). Due to the 
variation in rock types and rock magnetic properties a 
blanket treatment could not be prescribed for all sites. The 
minimum dispersion criterion (Irving et al, 1961) was 


applied to identify the best estimate of the mean remanence 
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direction at each site. Means were calculated ina 
hierarchial fashion. The values of D and I corresponding to 
the lowest value of a,, (highest Fisher k) were selected 
from among the sets of meaSurements at each site. In no case 
were fewer than three specimens from different cores finally 
selected to yield a D and 1, and in most instances all cores 
could be used with no rejections. At nine sites one or two 
prec anene exhibited anomalous behavior, either by 
Systematically approaching the group formed by companion 
cores but never quite reaching this goal, or by stubbornly 
retaining a divergent direction. These specimens were 
excluded from further analysis, but it is noteworthy that 
they represent only 5% of the total data set. Although in 
part subjective, it is felt that this procedure yields the 
best available estimate of the true mean paleomagnetic 
vector at each site. A series of field and laboratory notes 
have been prepared site by site in descending stratigraphic 
order as a reference for this collection of data and to 
clarify any peculiarities at a given site. These notes are 
presented in Appendices 1 and 2. 

The directional variance within each unit is generally 
small, with Fisher's precision parameter (Fisher, 1953) 
kedOO0Rin “S4eg0utiot S0@tiows ,Candekn<50 annoniynJecases. 
Another convenient measure of the quality of the directional 
data is provided by the angular shift occurring at the 
demagnetization step used in the final analysis. When sites 


with poor internal grouping are excluded, as explained 
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below, these angular shifts are quite small. Sixty-seven 
percent of specimens undergo a shift of 3° or less, the 
corresponding values for <5° and <10° being 79% and 92% 
respectively. Thus a small fraction of the specimens have no 
direct evidence of end-points, and their inclusion rests on 
agreement with cores from the same Site, as reflected by the 


appropriate site statistics. 
MAGNETOSTRATIGRAPHY FOR LEVEL MOUNTAIN 


The Stratocone Section, (PB)-Meszah Peak 

The essential magnetostratigraphic information for 
Level Mountain 1S Summarized in Table 11-1. The variation of 
VGP latitude as a function of elevation is presented as a 
magnetogram in Figure 11-2. 

Tt is evident that the stratocone section (PB) spans at 
least three polarity zones. In order of increasing age these 
are normal, reversed, normal. The uppermost normal 
magnetozone, Stratigraphically map unit 9, comprising 
hawaiite flows at the summit of Meszah Peak, can confidently 
be regarded as sampling the Brunhes (Epoch 1). These flows 
are the most recent on Level Mountain and field relations 
indicate that they postdate Quaternary continental 
glaciation (according to field mapping performed by the 
author and observations of Ostensoe, 1960). 

Sites PBP down through PBL are clearly reversed. The 


most straightforward correlation with the established 
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evel Mountain Paleomagnetic Remanence Data by Sites 
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elevation above sea level (feet) 

mean density in g/cc for cores from site 

% of NRM intensity remaining upon reaching endpoint 
field strength of last demag step upon reaching endpoint 
number of cores accepted (rejected) 

vector resultant of n unit vectors 

Fisher's precision parameter 

half-angle of cone of 95% confidence 

declination 

inclination 

latitude of equivalent VGP 

longitude of equivalent VGP 
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FAQGUBOW ila), 
Magnetogram for Level Mountain lavas showing polarity (VGP 
Lat) versus absolute stratigraphic elevation for the two 
Level Mountain sections. Lines are broken for erosional 
interval or other indicated time hiatus. 


MAGNETOGRAM FOR LEVEL MOUNTAIN 
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geomagnetic polarity time scale, Figure 11-3, is that these 
belong to the Matuyama (Epoch 2). It is uncertain whether 
this sequence of five reversed sites at Level Mountain lies 
above, below or misses the Jaramillo event entirely. It is 
not unreasonable to suppose that the short term normal 
events (Jaramillo, Gilsa, Olduvai) comprising less than 
0.3MY (McDougall and Chamalaun, 1966; Cox, 1969) could be 
missed by a random sampling of only a few flows spread over 
the length of the Matuyama (some 1.75MY). 

Sites PBY and PBK at the base of map unit 6 are normal, 
but of low precision, and have somewhat shallow 
inclinations, possibly representing a transition or 
excursion of the geomagnetic field. Flows at the base of map 
unit 6 have been dated (by fission track on pitchstones) at 
two locations on Level Mountain near the head of the 
Kakuchuya Valley. The dates. are 2.61 +0.12 and 2.46 +0.12MY 
(N. Briggs, pers. comm., 1978). These dates coincide within 
error to the Matuyama/Gauss boundary. It appears that the 
sequence of flows from PBL down through PBJ samples this 
transition, with sites PBY and PBK actually lying in the 
transition itself. 

Sites PBJ, PBI and PBH most likely lie in the upper 
portion of the Gauss (Epoch 3). Below these basalts is a 
Stratigraphic break and a lithological change to peralkaline 
trachytes. The next four flows have very shallow remanence 
directions, but the corresponding within- and between-site 


grouping is good. This group of flows (PBD, PBE, PBF, PBG) 
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FIGURE AT=S. 
Reference geomagnetic polarity time scale adapted from 
Mankinen and Dalrymple (1979) down to the base of the 
Thvera, and from MacDougall et al (1977) down into Epoch 
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represents either an excursion or a transition from an 
underlying normal to an unrecorded reversed polarity event 
lying between PBG and PBH. This could be the transition from 
the lower Gauss into the Mammoth event, the transition into 
the base of the Kaena, the transition out of the top of the 
Cochiti or the transition out of the top of the Nunivak. 
Without additional sampling or radiometric dating, this 
ambiguity cannot be resolved. The underlying three flows, 
PBC, PBB and PBA give normal directions, if somewhat 
shallow. Despite the shallow nature, the distribution of 
directions at PBA at least is very tight. Stratigraphically 
and petrochemically, site PBX belongs with map unit 5b and 
the sequence of basalts above, but is reversed. This site 
has no clear endpoint' but neither does it display random 
behavior to AF demagnetization. It is tentatively assigned 
to the top of the Gilbert (Epoch 4) implying that the base 
of map unit 5b samples the Gilbert to Gauss transition at 
3.40MY (Mankinen and Dalrymple, 1979). The reason that this 
second interpretation is not favored is that the overlying 
sequence of seven flows would all have to fit into the 
40,000 year normal interval between the Mammoth and the 


Kaena. These possibilities could be further assessed with 


' On AF demagnetization as the softer (and usually random) 
components of remanence are progressively removed, the 
remanence direction changes systematically approaching the 
true TRM. Once the TRM has been reached, further 
demagnetization produces no change in direction until all of 
the remanence has been removed. Sister specimens and other 
cores from the same site also approach this common TRM 
direction, which is termed the endpoint. 
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another stratigraphic section. There is unsampled 
Stratigraphy adjacent to PBX but it would have to be picked 
up laterally (westwards) along the northern Kakuchuya valley 
wall. The problem in the vicinity of the PB section is the 
morainal cover at this stratigraphic interval where the 
hanging valley enters the Kakuchuya. 

The underlying map unit 5a is the basal unit of the 
stratocone. Sites PBU and PBV (phonolites) and site PBW 
(comendite) are all normal. This fact, and the foregoing 
discussion, would place the three flows in one of the 
short-lived normal events of the Gilbert reversed epoch. 
Another pitchstone fission track date of 4.51 +0.26MY was 
obtained from this unit 5a further to the south across the 
Kakuchuya-Dudidontu ridge (N. Briggs, pers. comm.). Assuming 
this date to correspond to PBW, this site (within error) 
could correspond to the Sidufjalli normal event between 4.47 
and 4.32MY (Mankinen and Dalrymple, 1979), or possibly the 
underlying Thvera. One other possibility is that the 
sequence of sites iS not quite as old as unit 5a laterally, 
and that those flows of the PB section are younger than the 
date given, possibly placing them during either the Nunivak 
or the Cochiti event. 

In summary, for the stratocone section, stratigraphic 
and map unit divisions correspond fairly well to magnetic 
breaks. There are many possible correlations between the PB 
section and the geomagnetic polarity time scale. The favored 


interpretation is that this sequence of twenty-five lavas 
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Spans most of the recent four geomagnetic polarity epochs, 
commencing during the Gilbert at about 4.5MY and continuing 
up into the Brunhes, younger than 0.73MY. Alternatively, 
assuming that all of the normal and reversed polarity zones 
in the lower part of the stratocone have been sampled 
results in a minimum age eStimate for the base of the 
Stratocone between 3.40 and 3.15MY. In the favored 
interpretation, this sequence catches: the Gauss/Matuyama 
transition, the beginning of the transition into the Kaena 
or Mammoth event, and the transition from the Gilbert into 
the base of the Gauss. Eruptive activity appears to have 
been sparse and sporadic because the sampling indicates no 
evidence for the various events in the Brunhes and Matuyama 
Epochs. 

The Plateau Section (PA)-Little Tahltan Canyon 

The section of 25 flows sampled from the Little Tahltan 
Canyon crosses one of the best exposures of basalts which 
comprise the Level Mountain shield. This section is lower 
both in elevation and stratigraphic position than the 
Stratocone section described above. 

There is no radiometric data at all for the PA section 
and thus any correlation with the standard time-scale must 
be regarded as rather tentative. Furthermore, Mankinen and 
Dalrymple (1979) do not carry their polarity time scale 
beyond about 5MY. Taking the base of the PB section as 
representing the Sidufjall event, and simply matching the 


sequence of normal and reversed intervals determined by 
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McDougall et al (1977),see figure 11-3, leads to the 
conclusion that the lowest sites (PAX and PAY) fall near the 
middie of Epoch 6 (about 6MY). This is a minimum estimate 
for the onset of volcanism in the area. If "events" within 
the Gilbert epoch or Epochs 5 and 6 have been missed in this 
sampling, which seems probable, the base of the section will 
be correspondingly older. In this context it is reasonable 
to suppose that the eight flows (PAE to PAL) represent the 
lower part of Epoch 5 (5.50-5.83MY), not the very short 
normal interval at the top of Epoch 5 (5.34-5.37MY) which 
one-to-one matching implies. It also seems probable that the 
extremely short events at 6.12 and 6.39MY would not be 
recorded at Level Mountain. If this is so then sites PAQ to 
PAS may represent the normal interval 6.27 to 6.39MY, with 
Sites PAX and PAY consequently having an age of about 6.6MY. 
The frequency of volcanism during the plateau building 
is apparently greater than during the subsequent stratocone 
Stage. For the minimum estimate of the onset of volcanism in 
the area of about 6MY, the duration of plateau building is 
about 1.5MY and possibly as long as 2.1MY for the case of 
missed events. This translates to an average time interval 
between eruptions of sixty to eighty thousand years, or 
about one order of magnitude more frequent activity than 
during the stratocone stage. However both the field 
relations and the remanence vectors indicate that eruptive 
activity during plateau time was more clustered. From this 


viewpoint, the four plateau units may have been short lived 
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intense periods of activity separated by quiescent periods 
of about half a million years in duration. For this second 
interpretation the quiescent periods during the plateau and 
Stratocone stages are of about the same length, but the 
volumes of erupted material are more than an order of 
magnitude different. 

A Comparison of the Paleomagnetic Sections from Level 
Mountain and Mount Edziza 

Magnetograms for the three paleomagnetic sections from 
Mount Edziza (Souther and Symons, 1974) are shown in Figure 
11-4. Tentative correlations of Level Mountain and Mount 
Edziza eruptive sequences to the geomagnetic polarity-time 
scale are given in Figure 11-5. The directions for sites 
PBS, PBR and pnotaares with the uppermost flows of 
Stratigraphic sections A and C on Edziza. The uppermost 
flows at both centres are basaltic in nature. 

Sites PBP, PBO (map unit 8) and PBN (map unit 7) have 
steeper inclinations and more westerly declinations than 
most of the sites from map units 5, 7 and 9 of Edziza which 
are also interpreted to be of Matuyama age. While Souther 
and Symons (1974) have apparently caught the Jaramillo 
normal event at their site A6, no Similar intervening normal 
event has been sampled at Level Mountain. 

The directions for sites PBP, PBO and PBN are 
tentatively likened to the direction at Edziza's site A5 
which is placed below the Jaramillo, while the intermediate 


chemistry of flows PBP and PBO (benmoreites) is similar to 
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Figure 11-4, 
Magnetograms for the three Edziza Stratigraphic, sections. 
Data from Souther and Symons (1974). Elevation referenced to 
base of each section. Bar scale indicates map unit 
designation. Connecting lines are broken for time hiatus (Or 
Sampling gap) as inferred from Stratigraphic notes. 
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Figure 11-5, 
Suggested correlation for volcanism at Level Mountain and 
Mount Edziza with the reference polarity time scale. 
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trachybasalts jofipunit 7) lin sections (AsandiGtat BRdziza: 

The eruptive activity of map unit 6 (Level Mountain) 
appears to have been contemporaneous with glacial activity, 
as it 1s underlain by and includes tills, and is overlain by 
fluvioglacial sands and gravels. While the overlying unit 
could be Pleistocene, the underlying till implies a 
glactation n> the region vatsor.priorsto t2e5MY< The 
combination of information on glaciation at Level Mountain 
and Mount Edziza indicates that 3 or more episodes of 
continental galciation occurred between the Pleistocene and 
about 5MY. 

The stratigraphic interval including map units 5 and 6 
at Level Mountain is lithologically most similar to the 
sequence from B16 to B32 (Edziza unit 2). These sections are 
both characterized by abundant peralkaline flows and tuffs. 
It appears that the bimodal volcanism and stratocone stage 
at Edziza may preceed its equivalent at Level Mountain by as 
much as 0.5MY, with such activity beginning at Edziza during 
the reversed magnetozone of Epoch 5 or at the base of the 
Gilbert. It is also interesting to note that salic volcanics 
are included in the Edziza plateau, while at Level Mountain 
they are restricted to the Stratocone cap. There is a 
general age, magnetic and lithologic similarity for the 
stratocone section at Level Mountain and the younger 
portions of the 3 sections on Edziza. 

The lower plateau section (B) at Edziza is bracketed by 


K-Ar dates as having been erupted between 6.1 and 4.9MY 


Hag 
| Tah ver * aaa 
‘ F oa rt he 
a aie A ; Dp a “i 
et: : ; 


ie he iy Hansa anol aa cy, aan satel 


ie 2 ta anna) a phi ona 9 es aha oud ae . 
seen Ma bw Sita: a atelpetiae: chapel: Wott, Oy aise 
7 Hea \crel: + rhabe! ‘gt ‘ake th niad abi ab a 


- io 


ie vi wi cere yt nee i 2 Luh | pak ab rred anateen Fe 


my 
ds aoc: \ Des Oo eenk ag mca seg sd 


MES PAG) Dade ee ews dyad vot ante 


os. 


ral fy OS OTME DALY Bee. Tee del > ime, de pobist O° 
- : i” o 
eC The pris 46 ‘d raid, nt nolart esinbe . 
‘ays if nie 440 t tatu i) notes & etatae a 
5 ‘ 3 ) rm it 


ai) ‘Wa Sek- Davee ia me garg see eet, 


bie: a ok bias i jarboaltnis ei dbagnvaw | wl 
oy) Moet 4 (9 Sie ae ORR oF ate gent eons 
ret S ant cpm meri ee: spi: ee hawt nad oath b> a 
1 PALS fn aba & ¢ he hee Tebow is. aul inet? a4eaganly : 
4 bated 26 sitacaytape ‘a4 ean ynd shi 
3 dein Hird peteian srg ie NENG ae. iz a 
ae ints od de eee cova te apypezeopeay 6 sve? “a 
nBsniov. Shise tans yon. oo > alana ite ab ae bi +1 alates ant | 
, ingdet Leva ccstimih | dara. beasts son's at Brakes fated oie, ; 


Ae 
Wy See raw. ‘hea anal: oa eed ee 233 ots ‘gait? a 


ren ried a SS > palate ate Alas yao6 L S188ge 
FARMNe: BT DNs ieecleM, Lewes hes ned 3 3 wnosctas4s 
it at ee iinea't ast? ng. ania ttaoq. 
we hes odse: 4 26 5 Ee i, oa abe usgstodag naeot ont 
| nine peutiniad: a ast68 2A 


( 
aig yas 


s CM: Swe Cot eaweee 
f 


432 


(Souther and Symons, 1974). The best estimate from combined 
paleomagnetic and radiometric information is that the Edziza 
plateau building commences sometime during the lowest normal 
event of Epoch 6 sampling all three polarity zones of that 
epoch and the overlying reversed magnetozone of Epoch 5. The 
low latitude VGP site (PAW) could possibly represent a 
transition to the short normal zone at the top of Epoch 7. 
This indicates the possibility that the lowest sampled 
Stratigraphy at Level Mountain is older by at least 0.4MY 


than the age cited for the base of the Edziza plateau. 
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CHAPTER 12. SUMMARY AND CONCLUSIONS REGARDING THE 
PETROGENESIS OF THE LEVEL MOUNTAIN VOLCANICS AND THEIR 
RELATIONSHIP TO THE TECTONICS OF THE INTERMONTANE BELT. 

The most salient geological feature of the Stikine 
region of northwestern British Columbia is the Late Cenozoic 
volcanic activity. Since upper Miocene time more than 
1700km* of lavas have been erupted both as isolated flows 
and major volcanic edifices covering about 12% of the land 
area in the Stikine Volcanic Belt (Souther,1977b). Within 
this region are eleven major volcanic centres of which Level 
Mountain and Mount Edziza are the largest. Most of the 
Stikine lavas belong to the sodic alkali basalt series and 
are of continental affinity. At level Mountain there are two 
volumetrically dominant lava types. Alkali basalt flows and 
related basic lavas (sub-alkaline basalt, hawaiite, 
ankaramite) comprise an oval shield approximately 50km in 
diameter with an average thickness of 600m. The 
shield-building stage lasted 2MY or more, ending during the 
Pliocene at about 4.5MY ago. At that time a fundamental 
change occurred both in the chemistry of the erupted lavas 
and their physiographic expression. This second stage of 
volcanism at level Mountain, with much greater variability 
in the eruptive style (including flows, pyroclastics, domes, 
cupolas, dykes and hypabyssal stocks), resulted in the 
construction of a complex centrally-placed stratocone cap 
approximately 25km in diameter with a thickness of about 


250m, tWhilettheteneitertclassie sodicealkalicbasalt 
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differentiation series is present, the dominant lavas of the 
stratocone stage were peralkaline trachytes. Other 
compositions in Bee a cane order of abundance are: 
comendites, rhyolites, hawaiites, benmoreites, phonolites 
and mugearites. 

A variety of models for magma transport were calculated 
uSing one-atmosphere viscosity and liquidus determinations 
for the range of lava compositions. Basic magmas had very 
rapid ascent velocities in the range of 1 to 10cm/sec with 
salic magmas being about two orders of magnitude slower. The 
volumes of typical flows and other single eruption features 
weremwon the orderfoiri0e.stoet0 *kmos'Thatythe: Level 
Mountain range iS comprised of many such low volume and low 
viscosity flows argues for low overall rates of volcanicity 
and the relatively efficient transport of small aliquots of 
magma directly from the upper mantle. 

The isotopic and chemical evidence for the origin of 
the Level Mountain volcanics is straightforward. Basalts and 
inclusionsmwith¢primitavels’Suy*«Srevaluest (0. 7025ato 
0.70440) and primitive Pb isotope ratios (in the range of 
typical values for MORB and oceanic islands) argue for an 
upper mantle derivation from a heterogeneous source. Oxygen 
isotope values for basalts and inclusions 
(5.5 to 5.9°/,,.0'*O0 SMOW) confirm that the upper mantle 
beneath Level Mountain is normal with respect to oxygen. The 
shift to lighter isotopic compositions for the stratocone 


basalts, with respect to those on the shield stage, 
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indicates some heterogeneity for oxygen in the source region 
as well, and may indicate a shift in the depth of the 
melting zone as volcanism proceeded. Considering all of the 
Late Cenozoic lavas from the IMB, scatter in the O and Sr 
isotope Signatures, along with the variation in major and 
trace element ratios for lavas of similar evolution, argues 
for a very local scale heterogeneity in their upper mantle 
source regions. Somewhat lower abundances of K, Rb, U, Th 
and light rare earths for lavas of the IMB, compared to 
those from the source regions of alkalic magmas in general 
(White et al, 1979), may indicate that the upper mantle 
beneath the Cordillera has been depleted by prior melting 
events. Such depletion might have occurred through volcanism 
and plutonism accompanying the Mesozoic emplacement of 
allochthonous Cordilleran terranes (Monger and Price, 1979). 

The major and trace element chemical variations for all 
of the Level Mountain lavas are inconsistent with the 
derivation of the volumetrically important peralkaline lavas 
of the Stratocone stage via the classical fractional 
crystallization processes from an alkali basalt primary 
magma. All of the chemical variation is typified by two 
clusters which are interpreted to indicate that there are 
two principle magma types, basalt and peralkaline trachyte. 
The peralkaline salic magmas have some equally primitive 
lead, strontium and oxygen isotope values to the basalts, 
but the range in oxygen and strontium extends to higher 


values indicating a combination of secondary processes 
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including crustal contamination and selective open system 
behavior. The peralkaline salic magmas are traditionally 
interpreted to be derived by substantial fractional 
crystallization of a wollastonite-normative intermediate or 
basaltic parent (the orthoclase effect of Bailey and 
Schairer, 1966). At Level Mountain the major element 
variation of these peralkaline lavas (principally Na and K) 
is inconsistent with a derivation via any extensive alkali 
feldspar fractionation. This interpretation is also born out 
in the normalized rare earth patterns which are totally 
lacking in europium anomalies. The distribution of salic 
lavas about the thermal minimum in the liquidus surface 
indicates that they may be a related series of partial 
melts. Geobarometry estimates for trachytes uSing silica 
activitiesS give equilibration depths below 20km, in the 
lower crust. It is uncertain whether this reflects a depth 
of origin or rest point on the ascent path for the salic 
magmas. The high one-atmosphere temperatures determined in 
this study and the expected shape of the liquidus with 
increasing pressure (Bailey et al, 1974) suggest that the 
melting probably does not occur within the crust. 

The source compositions and P, T conditions for the 
Origin of the basaltic and peralkaline magmas are not very 
well constrained. Geothermometry and geobarometry 
calculations on lherzolite inclusions from the Stikine 
Volcanic Belt indicate equilibration depths in the range 


from 30 to 70 kilometers. Geobarometry estimates from silica 
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and alumina activities on Level Mountain lavas gives a 
Similar depth range. This agrees with data for central B. C. 
reported by Fiesinger and Nicholls (1977) and Fujii et al, 
1981. According to the recent synthesis of P and S wave 
residuals for the Cordillera, the upper mantle lacks a 
seismically resolvable high velocity lid, and these 
aforementionned depth estimates correspond to a range from 
the base of the crust down through the upper half of the low 
velocity zone (Wickens and Buchbinder, 1980). While the 
geobarometry estimates place the likely source region in the 
Stability field for spinel lherzolite, the lherzolite 
compositions from the IMB are too geochemically depleted 
(Tatsumoto, 1978) to give rise to extractable quantities of 
basaltic partial melts. Matrix calculations give maximum 
partial melts of 4% but less than 1% is far more common. 
Similar calculations using Level Mountain basalts and a 
pyrolite source give partial melts in the range from 5% to 
10% which is Higimendugh for melt extraction. The salic 
lavas require either a more enriched source or a different 
derivation mechanism. Peralkaline partial melts can only be 
derived from peralkaline bulk compositions, which is an 
unlikely situation in the crust or the mantle. One 
possibility is that the upper mantle or portions of it 
become metasomatized by such processes as gaseous transfer 
of incompatible elements from deeper in the mantle or by a 
series of partial melting episodes with inefficient melt 


extraction (Boettcher and O'Neil, 1980). This differentiated 
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upper mantle could contain pods of compositions approaching 
granulite and gabbro, which is what the partial melting 
calculations for deriving salic melts require. 

The cause of the Level Mountain and related Late 
Cenozoic volcanics of the IMB is an open question. Alkaline 
magmas are associated with extensional tectonics such as 
continental rift and oceanic fracture zones but they also 
occur on oceanic islands, seamount chains and in continental 
Margin environments (MacDonald, 1974; Bailey, 1977). Some of 
the structural characteristics of the IMB, as deduced from 
regional geophysics, resemble a continental rift setting. 
This structural similarity and the presence of alkaline 
volcanism, led Souther (1977b) to hypothesize an incipient 
continental rift. A drawback to this model is the lack of 
active seismicity in the IMB (Milne et al, 1978) and that 
the magnitude and rates of volcanism for the SVB are lower 
than Bast Africa by two orders of magnitude. The related 
alkaline volcanism of the Anahim Volcanic Belt led Bevier et 
al (1979) to point out a possible association with a tear 
fault in the subducted Juan de Fuca plate. There has been no 
subduction activity ingthe vicinity of the northern Canadian 
Pacific margin since approximately 20MY ago (Grow and 
Atwater, 1970). This Jack of a clear cut tectonic 
association led to a critical assessment the energy 
requirements for Level Mountain and the Stikine. 
Calculations on the amount of energy dissipated by volcanism 


indicate that thesbate (Cenozoic. volcanismrof the SVB could 
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be initiated and maintained by small changes (15%) in the 
total heat flux for the region. Changes of this magnitude 
could arise from the small tectonic stresses associated with 
the modern shear geometry of the adjacent continental margin 
(slip on the Queen Charlotte fault is about 6cm/yr) or from 
the re-equilibration of the thermal structure for the 
underlying upper mantle, with uplift and rebound of the 
geotherm following the completion of the subduction for the 


Kula plate in Mid-Tertiary time. 
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APPENDIX 1. FIELD AND REMANENCE NOTES ON LEVEL MOUNTAIN 


STRATOCONE SECTION 


COMMENTS (PB section decreases in elevation to 
south) 


PBT Hawaiite summit of Meszah Peak, map unit 9. Magnetic 


PBS 


PBR 


PBO 


PEE 


declination very erratic in site vicinity, possibility 
of lightening strikes, directions approach selected at 
different rates from different directions, inhomogeneous 
magnetization between cores, sister specimens often 
disagree until cleaned to >800 oersteds. The chosen 
direction is not a stable endpoint for all cores. The 
rejected ones were still coming in. Cores 2 through 5 
topographically oriented on Nuttlitude cone. 


Hawaiite next to highest flow on west shoulder of Meszah 
Peak between crater and dyke, map unit 9. While three 
cores approach a common endpoint, cores 4 and 5 do not, 
indicating the site to be inhomogeneously magnetized. 


Hawaiite highest prominent cliff outcrop on West 
shoulder of Meszah Peak as seen from south side, map 
unit 9. Abundant xenoliths, all cores approach a common 
endpoint in 200 to 400 oersted range. Restricting choice 
to cores 1,2 and 5 at 400 oersteds would improve 
Statistics without changing direction. 


Hawaiite cliff former below PBR map unit 9, abundant 
xenoliths, all cores approach common endpoint at 200 to 
400 oersteds. Restricting choice to cores 1,2 and 5 
improves statistics without changing direction. Below 
this elevation steep slopes, talus and scoria limit 
available outcrop. Additional section between Q and P 
could be collected on East face of Meszah Peak with aid 
of rope work. 


Benmoreite map unit 8. Narrow outcrop <3m wide in 
creek-snow chute on South face of Meszah Peak. Flows at 
this level are all thin and aside from a few blocky 
outcrops like this are scoriacequs, many thin trachytes 
interspersed, frequency decreasing upsection. All cores 
approach a common endpoint and stay there between 25 and 
400 oersteds. 
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Appendix 1. continued. 


PBO 


PBN 


PBM 


PBL 


PBY 


PBK 


PBJ 


PBH 


Benmoreite map unit 8, lowest blocky outcrop in snow 
chute, abundant mafic inclusions, below this point talus 
and scree cover resumes, all cores are stable at a 
common endpoint from NRM to 800 oersteds. Restricting 
choice to cores 1,2 and 3 improves statistics without 
changing direction. 


Trachyte lava tube low on south face of Meszah Peak, map 
unit 7A. All specimens are stable and coincident from 
NRM to 400 oersteds. 


Trachyte ledge, unit 7A at break in slope between lava 
tubes and tarn, all cores coincident from NRM to 400 
oersteds. 


Trachyte ledge unit 7A. Site inhomogeneously magnetized 
probably due to low and variable magnetic mineral 
content. Cores come in from different directions at 
different rates. Endpoint somewhat scattered. statistics 
can be improved without affecting direction by 
restricting choice to cores 1,2,3 and 4 at 400 oersteds. 


Trachyte flow massive columnar jointing, flow 46m thick, 
unit 7A, cores sampled from separate columns 2/3 of the 
way up the flow, site is transitional and scattered. No 
obvious stable endpoint is achieved. Intensity is low 
and remnance is soft. Statistics can't be much improved 
unless sample is dropped to 2 cores. 


Basal trachyte flow of unit 7A, near base of cliff 
outcrop formed by PBY. Laterally unit 6 is overlain by 
fluvioglacial/volcanofluvioglacial sands and gravels. 
Site 1S scattered and transitional. Some of scatter 
could be due to poor coring as all specimens were either 
in three pieces, short, or hole oriented. 


Basalt, unit 6B. Columnar jointed outcrop forms 
waterfall, good agreement of all specimens from NRM to 
400 oersteds. 


Hawalite, unit 6B. Weathering light grey very abundant 
phenocrysts of plagioclase and black augite, good 
agreement of all specimens from NRM to 400 oersteds. 


Hawaliite, unit 6B, with very abundant plagioclase and 
black augite. Direction chosen as stable endpoint, 
Statistics can be markedly improved without changing 
direction by eliminating core 3, which comes in more 
Slowly than the others and shows some disagreement 
between sister specimens, possible hematite and CRM 
component in core 3. 
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Appendix 1. continued. 


PBG 


PBF 


PBE 


PBD 


PBC 


FBG 


PBA 


PBX 


PBU 


Peralkaline trachyte, unit 6A. Uppermost of three 
Similar flows, all cores agree with good statistics but 
inclination is very shallow, transitional site. 


Peralkaline trachyte, unit 6A. Forms ledge and holds up 
flat profile in stream valley, phenocrysts of 
anorthoclase,fayalite and green pyroxene, all specimens 
agree from NRM to 800 oersteds with shallow inclination, 
transitional direction. 


Peralkaline trachyte flow, unit 6A. Holds up prominent 
waterfall, peculiar exfoliation textures, good agreement 
of all specimens but shallow inclination, transitional 
direction. From this point down only the prominent 
resistant flows are sampled, tuffs and agglomerates are 
omitted. Laterally in several places a till is present 
at this stratigraphic level. 


Vesicular browm weathered hawaliite, unit 5B, with 
plagioclase and augite phenocrysts. All specimens 
coincide at a shallow inclination transitional 
direction. These cores all contain appreciable maghemite 
both as phenocryst alteration and as veins and 
amygdules. 


Amygdaloidal alkali basalt flow, unit 5B. All cores 
achieve stable endpoints but site is somewhat scattered. 
Statistics can be improved without affecting direction 
by eliminating core 3 which was hole oriented. 


Alkali basalt with plagioclas and black pyroxene 
phenocrysts, unit 5B. Site somewhat scattered. Core 2 
was eliminated due to inhomogeneity inferred from 
disagreement between sister specimens. 


Hawaiite flow, unit 5B, stable endpoint for four 
Specimens 


Hawaliite rubbly outcrop west of stream cut, unit 5B. No 
exposure in creek at this point due to lateral moraine 
of Kakuchuya valley, site is definitely reversed. 
Eliminating core 2 improves statistics but does not 
affect direction. 


Phonolite flows cooled as a single unit. Slight primary 
CPyrdip cows £OuUN230-8, unit “5A. site, is somewhat 
scattered. 
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Appendix 1. continued 


PBV Phonolite flows, unit 5A. Weak site, all specimens but 
one have common NRM direction but scatter away randomly 
by 100 oersteds indicating a soft remnance. 


PBW Crystalline comendite flow, unit 5A. Stable direction 
for all cores from NRM to 100 oersteds. 
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APPENDIX 2. FIELD AND REMANENCE NOTES ON LEVEL MOUNTAIN 
PLATEAU SECTION 


SITE COMMENTS (Little Tahltan section decreases in 
elevation to south along E facing cliff) 


PAA Alkali basalt flow, map unit 4. Columnar outcrop, local 
topographic high back from cliff edge, possible 
lightening strikes, inhomogeneous site all cores 
approach the chosen endpoint but at different rates. 
Cores 3 and 5 were rejected as they were still coming in 
at 400 oersteds but became random by 600. 


PAB Alkali basalt flow, map unit 4, 3m thick columnar flow 
sampled im up from base. This was the next prominent 
outcrop 150m south of PAA, very coherent site, stable 
from NRM to 400 oersteds. 


PAC Hawaiite flow, map unit 4, abundant phenocrysts of 
plagioclase and black clinopyroxene, good site, all 
cores agree from NRM to 200 oersteds. 


PAD Alkali basalt, map unit 4. Site 1s reversed but 
Scattered, while individual cores attain stable 
endpoints (1 and 4), there is no common direction. Each 
core came from a different column and the flow was 
sampled over a distance of 15m laterally but this is not 
atypical. 


PAE Hawalite map unit 4, very tight site. All cores agree 
from NRM to 400 oersteds. 


PAF Hawaiite, map unit 4. No safe acess at this point in 
section so PAF was sampled 300m line of site to the 
south, good endpoint for all cores 100 to 200 oersteds. 


PAG Hawaliite, map unit 4. Good endpoint. All specimens 
Stable and coincident from 50 to 400 oersteds. 


PAH Hawaliite basal columnar flow of map unit 4 rests on 
yellowish green tuff horizon that delimits units 3 and 4 
over much of the southern plateau, stable direction, ail 
cores agree from NRM to 200 oersteds. 


PAI Basalt highest major cooling unit of continous exposure 
and thickness of unit 3, Spheriodially weathering coarse 
ophitic texture. The five cores here were sampled from 
only two columns, good site, all cores agree in 100 to 
200 oersted range. 
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APPENDIX 2. continued. 


PAJ 


PAK 


PAL 


PAM 


PAN 


PAP 


PAQ 


Alkali basalt, map unit 3, good site, all specimens 
agree in 100 to 200 oersted range. 


Alkali basalt, map unit 3, calcite and zeolites present, 
rejecting core 5 improves statistics without 
Significantly changing direction. 


Alkali basalt, basal flow of 3 that comprises a single 
columnar cooling unit, forms promontory north of stream 
notch in cliff. While the overlying section was sampled 
along 2km of cliff margin, the underlying section 
follows the stream notch in this cliff and is nearly 
vertical. All cores for this suite attain a good common 
endpoint between 100 and 200 oersteds. 


Basalt flow with thompsonite and chabazite in fractures. 
This is third flow down from PAL due to poor acess, map 
unit 3. The intervening flows may contain the transition 
from the reverse polarity of PAL. The stable endpoint is 
achieved by all cores at 200 oersteds, the direction is 
a little shallow making this a transitional site. 


Alkali basalt unit 3, lower flow of pair that comprises 
a Single columnar cooling unit immediately below PAM, 
Site somewhat scattered similar shallow reversed 
direction to PAM, statistics can be improved by 
restricting choice to cores 2,3 and 5 but direction does 
not change. 


BaSalt unit 3, xenocrysts of forsterite and brown 
spinel? ;Mqoodesitentighes distribution for 125 .to-£00 
oersteds on all but core 4, which is consequently 
rejected. 


Basalt columnar outcrop, unit 3, overlies carbonized 
wood and peat, Site reversed and somewhat scattered best 
agreement in direction at 200 oersteds but core 4 still 
coming in. 


Alkali basalt, lowest flow in unit 3, upwards to PAP 
there are two flows and a rubbly Lahar which were not 
Sampled due to poor access on cliff face, this sampling 
gap coincides with a polarity change, this site is 
normal with good agreement for all cores above 200 
oersteds despite the fact that PAQ is mapped in unit 3 
and seems to belong there chemically as well, its 
magnetic orientation appears to go better with 
underlying PAR and PAS. The time interval between units 
2 and 3 is probably small. 
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APPENDIX 2. continued 


PAR 


PAS 


PAT 


PAU 


PAV 


PAW 


Alkali basalt separated from PAQ by an orange-buff tuff 
horizon mapped as the unit 2-3 boundary. There is good 
agreement for all cores by 200 oersteds except #2 which 
did not have a sun bearing and the two topographic 
Sitings disagreed. This core also did not achieve an 
endpoint so it should be rejected. 


Transitional hawaiite-basalt map unit 2, Two flows 
between here and PAR were omitted due to cover inside 
the notch and no access on the cliff face, core 5 had to 
be rejected because it disagreed with all other cores 
and was more intensely magnetized by 2 to 3 orders of 
magnitude, otherwise this is a good site. 


Vesicular altered hawaiite (?) flow map unit 2, two 
flows omitted between here and PAS. From here down cliff 
exposure required that safety ropes be used during 
drilling and orientation. This site is quite scattered, 
all cores achieve endpoints but only 1 and 2 coincide, 
the site is reversed but probably has a later CRM 
component related to maghematite alteration as seen 
petrographically, and in failure of core to saturate 
even by 5 kilogauss in IRM test. 


Basalt, map unit 2, massive columnar flow >5m thick, one 
flow between here and PAT not sampled, most cores 
coincide from NRM to 400 oersteds. 


Basalt lowest massive flow in map unit 2. This columnar 
cooling unit 1S comprised of four individual flows, one 
cooling unit between here and PAU not sampled. Cores 2 
and 3 were drilled through fractures and became broken 
in several pieces and had to be hole oriented. They 
Scatter somewhat from other 3 cores which show good 
agreement by 200 oersteds. 


Alkali basalt highest cooling unit in map unit 1, the 1 
to 2 boundary, mapped at the top of PAW because of 
distinctive colour and massive appearance of underlying 
flows. However there is no tuff horizon, paleosol or 
other obvious time break here. This site has an unusual 
Shallow inclination but is normal compared to adjacent 
flows above and below. It is also low intensity for 
basalts, probably represents an excursion Since long 
time break is not inferred. 
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APPENDIX 2. continued. 


PAX Transitional basalt massive cooling unit typical of map 


PAY 


unit 1, good agreement for all cores from NRM to 200 
oersteds. 


Basalt massive basal flow of map unit 1. Below this 
level there is no further outcrop nearby, only several 
hundred feet of rubbly cold agglomerates 
(volcano-fluviatile). The Little Tahltan River here 
appears to be resequent to a Tertiary stream valley of 
moderate proportions. Bedrock 1S porphyry similar to 
outcrop on Kaketsa Mtn. and of the Sloko volcanics. 
Although all cores agree with high Fisher k the best 
direction is probably given by cores 1 and 4 at 100 to 
400 oersteds as other cores approach and pass through 
this region. 
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